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The  atomic  structures  of  four  Ca40+xMg25Cu35  x  (X  =  0,  5,  10,  and  20  at.  %)  ternary  metallic 
glasses  have  been  determined  using  a  synergistic  combination  of  neutron  diffraction,  ab  initio 
molecular  dynamics  (MD)  simulation,  and  constrained  reverse  Monte  Carlo  modeling.  It  is 
described  as  close-packing  of  efficiently  packed  Cu-centered  clusters  that  have  Ca,  Mg,  and  Cu 
atoms  in  the  first  coordination  shell.  The  close-packed  arrangement  of  the  clusters  provides  a 
characteristic  medium  range  order  in  these  alloys.  An  average  coordination  number  (CN)  of  10 
(with  about  5  7  Ca,  2  3  Mg,  and  1  2  Cu  atoms)  is  most  common  for  the  Cu-centered  clusters.  The 
average  coordination  numbers  around  Mg  and  Ca  are  12  13  (~6  8  Ca,  3  Mg,  and  1  4  Cu)  and 
13  15  (7  9  Ca,  3  4  Mg,  and  2  5  Cu),  respectively,  and  they  are  composition  dependent.  Strong 
interaction  of  Cu  with  Mg  and  Ca  results  in  pair  bond  shortening.  Icosahedral  short  range  order 
does  not  dominate  in  these  amorphous  alloys,  although  polytetrahedral  packing  and  five-fold  bond 
configurations  resulting  in  pentagonal  bi-pyramids  have  been  found  to  be  the  most  common  nearest 
atom  configurations.  ©2012  American  Institute  of  Physics.  [http://dx.doi.org/10.1063/L4729450] 


I.  INTRODUCTION 

It  is  well  established  that  Ca-Mg-based  metallic  glasses 
have  very  good  glass  forming  ability  (GFA). 1  6  Based  on 
two  simple  metal  elements,  these  glasses  have  properties  that 
distinguish  them  from  transition-metal-based  bulk  metallic 
glasses  (BMGs),  such  as  very  low  Young’s  and  shear  moduli 
that  are  comparable  with  the  moduli  of  human  bones,7  8  very 
low  density  (1. 6-2.8  g/cm3),9  and  strong  relaxation  dynamics 
of  the  super-cooled  liquid.10  Glass  forming  ability,  as  well  as 
physical  and  mechanical  properties  of  metallic  glasses,  is 
believed  to  depend  on  the  type  of  short  range  order  (SRO) 
and  medium  range  order  (MRO)  of  alloying  elements  in  the 
amorphous  structure.  It  is,  therefore,  interesting  to  know  how 
the  amorphous  structure  of  Ca-Mg-based  BMGs  differs  from 
transition-metal-based  BMG  structures. 

Structural  analysis  of  metal-metal  BMGs  has  mainly 
focused  on  transition  metals  glasses  such  as  Zr-Cu,11  11 
Zr-Pt,15  and  Zr-Cu-AI. 16,17  The  results  indicate  that  the  pack¬ 
ing  of  atoms  in  these  materials  is  not  random  but  is  strongly 
influenced  by  chemical  interactions.  In  particular,  icosahe¬ 
dral  SRO  has  been  identified  and  its  importance  on  relaxa¬ 
tion  dynamics,  GFA,  and  mechanical  properties  of  these 
BMGs  has  been  discussed.  On  the  other  hand,  transition 
metal  metalloid  glasses  such  as  Fe-C  and  Ni-B, 18-20  and 
many  A1  transition  metal  glasses  with  high  A1  content,21  are 
not  characterized  by  icosahedral  SRO.  These  observations 
indicate  that  metallic  glasses  can  have  various  types  of  SRO.-- 


a,Author  to  whom  correspondence  should  be  addressed.  Electronic  mail: 
oleg.senkov@wpafb.af.mil.  Telephone:  937  255  4064. 


The  atomic  structure  of  Ca60MgxZn4o  x  BMGs,  identi¬ 
fied  with  the  use  of  x-ray  and  neutron  diffraction  and  reverse 
Monte  Carlo  (RMC)  simulation,23  shows  no  icosahedral 
SRO.  The  basic  building  blocks  of  the  amorphous  structure 
in  these  glasses  are  Mg-  and  Zn-centered  clusters  with  pri¬ 
marily  Ca  atoms  in  the  first  coordination  shell.  The  average 
coordination  number  (CN)  around  Zn  is  9  10  [~6-7  Ca  +  3 
(Mg  +  Zn)],  with  dominant  (0,3, 6,0)  and  (0,2, 8,0) 
polytetrahedral-type  clusters,  and  the  CN  around  Mg  is 
11  12  [~7  8  Ca  +  4  (Mg  +  Zn)],  with  a  large  fraction  of 
(0,2,8, 1)  and  (0,2, 8, 2)  clusters.  Ca  atoms  have  an  average 
CN  of  13.  Ab  initio  MD  simulations  support  these  results. 
Although  full  icosahedra  do  not  control  the  atomic  structure 
in  Ca-Mg-Zn  BMGs,  five-fold  bonds  in  the  form  of  pentago¬ 
nal  bi-pyramids  are  the  most  populous  structural  units.2’  A 
large  fraction  of  five-fold  bonds  and  the  lack  of  icosahedral 
SRO  has  also  been  found  in  the  MD  simulated  amorphous 
structure  of  Mg-Cu  alloys.-4  The  average  CN  around  Cu  and 
Mg  is  composition  dependent,  changing  from  less  than  10 
with  no  icosahedral  clustering  to  above  12  with  increasing 
Cu  content.  Only  on  the  Cu-rich  side,  where  the  CN  is  close 
to  12,  does  a  predominance  of  five-fold  bonds  indicate  Cu- 
centered  icosaherdal  SRO,  which  is  similar  to  the  Zr-Cu 
case.11-14 

The  composition  dependence  of  GFA  in  Ca-Mg-Cu 
glasses  has  been  well  established.5'6  The  GFA  and  relaxation 
dynamics  of  some  of  these  alloys  are  similar  or  even  slightly 
better  than  those  of  the  best  Ca-Mg-Zn  BMGs.910  At  the 
same  time,  Ca-Mg-Cu  BMGs  have  better  corrosion  resist- 
ance  ’  and  seem  to  be  stronger  than  Ca-Mg-Zn  BMGs. 

In  addition,  strong  covalent-like  interactions  between  Cu-Mg 
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and  Cu-Ca  have  been  noticed.5  It  is,  therefore,  useful  to  com¬ 
pare  the  amorphous  structures  of  Ca-Mg-Cu  and  Ca-Mg-Zn 
BMGs.  The  ternary  Ca-Mg-Cu  system  has  clear  advantage 
compared  to  other  glass  forming  ternary  systems,  such  as 
Zr-Ti-Cu  or  Zr-Al-Cu,  because  the  very  different  atomic 
radii  of  Ca  (~1.9A),  Mg  (1.6  A),  and  Cu  (1.25  A)  make  the 
separation  of  the  partial  distribution  functions  relatively  easy 
and  unambiguous. 

In  the  present  work,  the  atomic  structures  of 
Ca4o+xMg25Cu35  x  BMGs  were  studied  using  a  synergistic 
combination  of  neutron  diffraction,  ab  initio  (quantum)  mo¬ 
lecular  dynamics  (QMD)  simulation,  and  constrained  reverse 
Monte  Carlo  (CRMC)  modeling.  First,  six  partial  radial  distri¬ 
bution  functions  (PRDFs)  were  determined  for  every  of  the 
studied  alloys  using  a  self-contained  QMD  simulation,  which 
requires  no  input  from  experiment.  The  simulated  PRDFs 
were  then  used  to  calculate  the  total  radial  distribution  func¬ 
tions  (RDFs)  and  compare  with  the  experimentally  measured 
RDFs.  Excellent  agreement  of  the  experimental  and  QMD- 
simulated  RDFs  demonstrated  that  the  QMD-simulated 
PRDFs  are  consistent  with  and  can  be  used  in  addition  to  the 
experimental  neutron  data  to  resolve  the  amorphous  structure. 
Finally,  the  3D  atomic  structures  of  the  studied  amorphous 
alloys,  which  are  consistent  with  the  experimental  diffraction 
data  and  QMD-simulated  PRDFs,  were  reconstructed  using 
CRMC  modeling.  The  simulated  atomic  structures  were  statis¬ 
tically  analyzed  using  pair  and  three-body  correlation  func¬ 
tions,  Voronoi  tessellation,  and  nearest  neighbor  approaches. 
These  results  were  compared  with  an  earlier  analysis  of  the 
atomic  structure  of  Ca-Mg-Zn  metallic  glasses.26 


II.  EXPERIMENTAL  PROCEDURES 

Fully  amorphous  samples  of  compositions,  Ca40+xMg25 
Cu35  X’  where  X=0,  5,  10,  and  20  at.  %,  were  prepared  by 
melt-spinning  in  the  form  of  ribbons  as  described  in  Ref.  23. 
The  densities,  pQ,  of  the  amorphous  alloys  were  measured 
with  a  helium  pycnometer  AccuPyc  1330  V1.03  and  the  val¬ 
ues  (in  g/cm3  and  atoms/A3)  are  given  in  Table  I.  Neutron  dif¬ 
fraction  experiments  were  conducted  at  room  temperature 
under  vacuum  using  the  general  materials  (GEM)  diffractome¬ 
ter  at  the  ISIS  high-intensity  pulsed  neutron  source  (Ruther¬ 
ford  Appleton  Laboratory,  Didcot,  UK).27  GEM  has  eight 
detector  banks  that  collect  data  over  a  wide  range  in  Q  (from 
0.1  to  100  A  *),  where  Q  —  An  sin0/A  is  the  magnitude  of  the 
scattering  vector  for  a  neutron  of  wavelength  A  scattered  at  an 
angle  20.  Just  prior  to  the  neutron  diffraction  experiment,  the 
ribbon  samples  were  crushed  into  powder  and  loaded  into 


TABLE  I.  Density  (in  g/cm3  and  atoms/A3)  of  Ca  Mg  Cu  amorphous  alloys 
studied  in  this  work. 


Alloy 

Density 

g/cm3 

Atoms/A3 

Ca40Mg25CU35 

2.936  ±  0.007 

0.03987 

Ca45Mg25Cll30 

2.673  ±  0.005 

0.03728 

Ca5oMg25Cll25 

2.439  ±  0.003 

0.03497 

Ca60Mg25Cui5 

2.039  ±  0.003 

0.03097 

10.3  mm  diameter  cylindrical  vanadium  containers.  The  thick¬ 
ness  of  the  container  walls  was  25  pm.  Data  reduction,  correc¬ 
tion,  and  analysis  were  carried  out  using  the  program  gudrun 
(Ref.  28)  and  the  atlas  suite  of  programs.29  This  yielded  the 
neutron  total-scattering  structure  factors  (SSF), 

3  i  3 

s(Q)  =  J2 cicMAMQ)  -  ']  =  -  J2 ctf-  (1) 

ij=  1  ‘=1 


Here,  is  the  differential  neutron  cross-section  per  unit 
solid  angle  Q  for  the  3-component  alloy,  c,  and  bt  are, 
respectively,  the  atomic  fraction  and  the  coherent  bound  neu¬ 
tron  scattering  length  of  element  i,  Ajj(Q)  are  the  Faber- 
Ziman  partial  structure  factors,30  and  ,  c,b2  is  the 

total-scattering  cross  section  of  the  alloy.  Each  S(Q )  was 
Fourier  transformed  to  give  a  real  space  neutron  total  RDF, 
G(r),  defined  as31 


G(r)  =  J2  c‘cMi[gij(r )  -  !] 
*J=  1 


oo 

1 

2n2rp0  . 

0 


QS(Q)L(Q)Sin(rQ)dQ. 


(2) 


Here,  pa  is  the  average  density  of  the  alloy  (in  atoms  per  A3), 
r  is  a  distance  from  an  average  origin  atom  in  the  amorphous 
structure,  and  gifjr)  are  the  PRDFs,  which  are  defined  as 


8ij(r)  = 


1  dnjj(r) 
Anr2Cjpa  dr 


(3) 


where  d/;^  are  the  number  of  elements  of  type  j  between  dis¬ 
tances  r  and  r  +  dr  from  an  element  of  type  i.  The  Lorch 
modification  function,  2  L(Q),  with  a  maximum  momentum 
transfer  Qnuix  of  25  A  1  was  used  to  reduce  termination  rip¬ 
ples  in  the  Fourier  transform  (Eq.  (2)).  Combining  Eqs.  (1) 
and  (2)  gives  the  following  relation  between  g,j(r)  and 
Aij{Q): 


2 n2rp0[gij{r)  -  1] 


OO 

Q[AU(Q)  -  mQ)MrQ)dQ-  (4) 

o 


Using  Eq.  (3),  the  partial  coordination  number  of  an  element 
j  in  the  first  coordination  shell  of  the  element  i  is  determined 


^max 

CNij  =  4ncjpa 

l  min 


gif2  dr. 


(5) 


Here,  rmin  and  rmax  are  the  positions  of  the  start  and  end  of 
the  first  peak  in  the  respective  gij(r). 

QMD  simulation  of  the  Ca-Mg-Cu  amorphous  structures 
was  conducted  using  the  Vienna  ab  initio  simulation  package 
(v asp). 33  The  projector  augmented-wave  (PAW)  method34'35 
and  Perdew-Burke-Ernzerhof  (PBE)  functional36  were  used 
to  describe  interacting  valence  electrons.  Cubic  simulation 
boxes  with  periodic  boundary  conditions  contained  200 
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atoms,  and  the  box  size  was  determined  by  the  density  of  the 
modeled  material.  The  liquid  was  equilibrated  at  1000  K, 
then  was  quenched  to  300  K  in  100  K  temperature  steps, 
equilibrating  at  each  temperature  for  3000  steps.  Another 
3000  consecutive  configurations  were  further  sampled  at  300 
K  with  a  time  step  of  5  fs,  and  an  averaged  PRDF  was  calcu¬ 
lated  from  these  configurations,  thus  taking  into  account 
thermal  vibrations  and  assuming  no  structural  relaxation 
within  15  ps.  The  QMD  model  provided  a  complete  set  of 
six  simulated  PRDFs,  gijQMD(r),  well  agreed  with  the  experi¬ 
mental  diffraction  data,  for  every  studied  alloy. 

To  expand  the  volume  of  the  simulated  structures  and  bet¬ 
ter  fit  to  the  experimental  diffraction  data,  alloy  density,  and 
closest  approach  constraints,  the  CRMC  simulation  technique 
described  in  detail  in  Ref.  37  was  used.  Simulation  boxes  with 
periodic  boundary  conditions  contained  19200  atoms,  and  the 
box  volume  was  determined  by  the  density  of  the  material. 
The  input  data  for  the  CRMC  simulation  were  S(Q),  G(r),  and 
six  g,yQMD(r).  Renonnalization  of  g/,QMD(V)  was  allowed  to 
correct  possible  QMD  simulation  errors.  More  than  10  million 
accepted  random  atom  moves  occurred  for  each  sample  during 
the  CRMC  simulation  before  the  final  configurations  provided 
acceptable  fits  to  the  experimental  and  QMD-simulated  data. 
This  assured  total  independence  of  the  modeled  structures  on 
the  initial  configurations  of  randomly  distributed  atoms. 

The  simulated  structures  were  statistically  analyzed 
using  pair  and  three-body  correlation  functions,  Voronoi  tes¬ 
sellation,  and  nearest  neighbor  approaches,  which  allowed 
calculation  of  local  structural  features  such  as  total  and  par¬ 
tial  CN,  triplet  angle  correlations,  type  and  distribution  of 
characteristic  coordination  polyhedra,  and  local  packing 
fraction  SRO  and  MRO. 

III.  RESULTS 

A.  Neutron  diffraction  analysis 

The  experimental  SSFs,  S(Q),  for  the  four  samples  are 
shown  in  Figure  1(a)  and  the  RDFs,  G(r),  are  shown  in 
Figure  2(a).  The  S(Q)  curves  damp  very  quickly  with  Q  and 
no  oscillations  are  apparent  beyond  ~15  A  ,  indicating  that 
the  glasses  have  a  wide  range  of  interatomic  distances.  As 
the  Cu  concentration  increases,  the  first  sharp  diffraction 
peak  (FSDP)  shifts  to  higher  Q  and  broadens  (Figure  1(b)). 
The  RDFs  show  that  an  increase  in  Cu  content  shortens  the 
average  interatomic  distance  and  decreases  the  distribution 
width  of  the  interatomic  distances  (Figure  2(b)).  Thus,  the 
shift  and  broadening  of  the  FSDP  with  an  increase  in  the  Cu 
content  can  be  interpreted  as  arising  from  the  shortening  of 
the  average  interatomic  distance  and  narrowing  of  the  distri¬ 
bution  of  interatomic  distances,  respectively.  There  is  also  a 
pre-peak  at  Qss  1.2  A  1  in  the  S(Q)  curves,  which  indicates 

1 5  38  39 

the  presence  of  medium  range  order  in  these  alloys.  ’  ’ 

Systematic  changes  of  the  shape  of  the  first  RDF  peak 
with  composition  are  observed  (Figure  2(a)).  In  particular, 
the  first  RDF  peak  in  Ca60Mg25Cu15  has  two  maxima,  one  at 
r  =  3. 1 1  A  and  another,  higher  intensity,  at  r  =  3.63  A.  As  the 
amount  of  Cu  increases,  the  relative  intensity  of  the  first 
maximum  increases  and  its  position  slips  to  lower  r-values 
(Figure  2(b)),  while  the  second  maximum  degenerates  into  a 
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FIG.  1.  (a)  Experimental  total  neutron  scattering  structure  factors  for 
Ca6o  xMg25Cui5+x  metallic  glasses.  The  curves  are  separated  by  shifting 
along  the  vertical  axis  for  better  view,  (b)  Dependence  on  the  concentration 
of  copper  of  the  positions  of  the  start  (1st  Minimum),  maximum  (Peak),  and 
end  (2nd  minimum)  of  the  first  peak,  as  well  as  the  peak  width. 
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small  bump  in  alloys  with  25%  and  30%  Cu  and  disappears 
in  the  alloy  with  35%  Cu.  At  the  same  time,  an  additional 
maximum  starts  to  develop  at  r  «  2.63  A,  preceding  the  first 
maximum,  in  the  alloys  with  30%  and  35%  Cu.  Such  a 
strong  concentration  dependence  of  the  shape  of  the  first 
RDF  peak  is  apparently  due  to  an  increase  in  the  fraction  of 
shorter  Cu-Cu  and  Cu-Mg  bond  pairs  and  a  decrease  in  the 
fraction  of  longer  Ca-Ca  and  Ca-Mg  bond  pairs  with  an 
increase  in  the  amount  of  Cu. 

B.  Ab  initio  molecular  dynamic  simulation 

Figure  3  shows  the  six  QMD-simulated  PRDFs  for  the 
Ca40Mg25Cu35,  Ca5oMg25Cu25,  and  Ca60Mg25Cu15  alloys. 
The  PRDFs  are  rather  smooth,  since  they  are  averaged  over 
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FIG.  2.  (a)  Total  radial  distribution  functions  for  Ca6o  xMg2sCui5+x  metal 
lie  glasses.  The  curves  are  separated  by  shifting  along  the  vertical  axis  for 
better  view,  (b)  Dependence  on  the  concentration  of  copper  of  the  positions 
of  the  first,  second,  and  third  peaks  and  the  width  of  the  first  and  second 
peaks. 
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3000  instantaneous  configurations.  The  effect  of  composition 
on  the  width  and  position  of  the  first  PRDF  peak  is  weak  rel¬ 
ative  to  the  effect  on  the  positions  and  intensities  of  the 
higher  order  peaks.  This  indicates,  to  the  first  order,  that 
bond  distances  between  respective  atom  pairs  in  the  first 
coordination  shell  are  not  affected  by  composition  and  that 
the  partial  coordination  number  of  element  j  around  element 
/  is  linearly  proportional  to  the  concentration  c,  and  alloy 
density  p0  (see  Eqs.  (3)  and  (5)).  From  the  QMD-simulated 
PRDFs  and  using  Eq.  (2),  simulated  total  RDFs  were 
obtained  for  every  alloy  and  compared  with  the  respective 
experimental  neutron  RDFs  (Figure  4).  A  very  good  match  is 
seen.  Since  QMD  simulations  were  done  completely  inde¬ 
pendently  from  the  neutron  diffraction  experiments,  these 


good  fits  indicate  that  the  simulated  PRDFs  are  consistent 
with  and  can  be  used  in  addition  to  the  experimental  neutron 
S(Q)  and  G(r)  to  resolve  the  amorphous  structures. 

The  average  partial  and  total  coordination  numbers 
around  Ca,  Mg,  and  Cu  atoms,  from  the  QMD  simulations, 
are  given  in  Table  II  and  the  most  frequent  (mode,  rjjode), 
weighted  average  (mean,  r^ean),  and  cut-off  (maximum, 
/[“'  °fi)  distances  between  atom  pairs  in  the  first  coordina¬ 
tion  shell  are  given  in  Tables  III  V,  respectively.  The 
weighted  average  distances,  r^ean,  between  the  /  and  j  atoms 
were  calculated  using  the  following  equation:-1 


Mean 

y 


rijgij{r)dr 


giM)dr. 


(6) 


C.  Constrained  reverse  Monte  Carlo  simulation 

The  QMD-simulated  PRDFs,  gyQMD(r),  were  used  as  soft 
constraints  in  CRMC  simulations  to  force  the  local  atomic 
order  to  be  consistent  with  experimental  S(Q)  and  G(r)  and 
simulated  g^QMD(r).  Renormalization  of  g^QMD(r)  during 
CRMC  allowed  minor  adjustment  of  the  QMD  simulation  to 
give  much  better  fits  of  the  CRMC-simulated  S(Q)  and  G(r)  to 
the  experimental  functions  (Figure  5).  The  CRMC  PRDFs  are 
shown  in  Figure  6.  Similar  to  QMD  results,  the  first  peak  posi¬ 
tion  almost  does  not  depend  on  alloy  composition,  however, 
the  peak  intensity  does.  Using  these  PRDFs  and  Eq.  (3),  the 
partial  coordination  numbers,  as  well  as  the  mode,  mean,  and 
maximum  cut-off  distances  between  the  pair  atoms  in  the  first 
coordination  shell,  were  calculated  and  the  results  are  given  in 
Tables  II  V,  respectively. 

1.  Voronoi  analysis  of  Ca-Mg-Cu  amorphous  structure 

Voronoi  analysis  was  used  to  study  local  atom  environ¬ 
ments  beyond  the  PRDFs  results.-1'40  In  Voronoi  analysis, 
atoms  that  have  common  faces  in  their  Voronoi  polyhedra 
are  nearest  neighbors  and  a  coordination  polyhedron  with 
vertices  at  the  first-shell  atom  positions  and  edges  coinciding 
with  the  interatomic  bonds  is  defined  for  any  chosen  atom 
i.4 1  Each  coordination  polyhedron,  also  called  an  /-centered 
cluster,40’42’43  is  associated  with  the  respective  /-centered 
Voronoi  polyhedron  and  can  be  assigned  a  Voronoi  signature 
( v3,v4,v5,v6 )•  F°r  a  Voronoi  polyhedron,  vm  is  the  number  of 
faces  containing  m  edges;  while  for  the  respective  coordina¬ 
tion  polyhedron,  vm  is  the  number  of  vertices  common  to  m 
polyhedron  edges  (or  faces).44  In  the  latter  case,  m  is  also 
called  the  vertex  coordination.  Those  /-centered  coordination 
polyhedra  with  the  same  Voronoi  signature  are  considered  to 
be  topologically  equivalent  (even  though  they  may  not  be 
identical)  because  they  can  be  transformed  into  each  other 
without  changing  the  number  of  vertices  and  connecting 
edges.  In  addition  to  assigning  topology,  the  Voronoi  sig¬ 
nature  also  defines  the  total  coordination  number  (C/V,)  of 
the  /-centered  cluster  as  CNj=^2vm.  Clusters  with  the 
same  Voronoi  signature  can  be  chemically  different. 
Therefore,  in  addition  to  the  Voronoi  signature,  partial 
coordination  numbers,  i.e.,  the  number  of  atoms  of 
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FIG.  3.  QMD  simulated  (a)  Cu  Cu,  (b) 
Mg  Cu,  (c)  Mg  Mg,  (d)  Ca  Cu,  (e)  Ca 
Mg,  and  (f)  Ca  Ca  PRDFs  for 
Ca4oMg25Cu35,  Ca5oMg25Cu25,  and 
Ca60Mg25Cui5  metallic  glasses. 


different  species,  should  also  be  known  for  a  more  com¬ 
plete  description  of  the  SRO  in  the  amorphous  structure. 
The  partial  and  total  coordination  numbers  obtained  with 
Voronoi  tessellation  of  the  CRMC-simulated  amorphous 
structure  of  Ca40+xMg25Cu3g  x  alloys  are  given  in  Table 
II.  They  are  in  good  agreement  with  those  obtained  with 
the  use  of  the  CRMC-simulated  PRDFs. 

Figure  7  illustrates  the  types  and  fractions  of  Ca-,  Mg-, 
and  Cu-centered  coordination  polyhedra  in  the  CRMC- 
simulated  amorphous  structures.  Although  many  types  of  coor¬ 
dination  polyhedra  are  present,  the  most  common  clusters  are 
(0,2, 8, 4)  for  Ca-centered,  (0,2, 8, 2)  for  Mg-centered,  and 
(0,3, 6,0)  and  (0,2,8, 1)  for  Cu-centered  clusters.  The  fractions 
and  distributions  of  the  clusters  depend  on  alloy  composition. 
Figure  8  shows  the  fractions  of  3-,  4-,  5-,  and  6-coordinated 
vertices  in  the  Ca-,  Mg-,  and  Cu-centered  clusters  in  the 
Ca-Mg-Cu  amorphous  alloys.  The  five-coordinated  vertices 
dominate  in  all  clusters,  which  is  an  indication  that  pentagonal 


bi-pyramids  are  the  main  building  blocks  in  these  amorphous 
alloys.  Taking  into  account  that  the  pentagonal  bi-pyramid 
consists  of  five  tetrahedra,  one  may  conclude  that  polytetrahe- 
dral  packing  plays  an  important  role  in  the  formation  of  the 
amorphous  structure  in  the  Ca-Mg-Cu  alloys. 

Radical  Voronoi  tessellation  gives  the  Voronoi  cell  vol¬ 
ume  (VCV)  of  each  individual  atom.  The  average  VCV  val¬ 
ues  of  Ca,  Mg,  and  Cu  atoms  are  given  in  Table  VI.  The 
VCV  is  higher  for  larger  atoms.  The  VCV  of  Ca  and  Cu 
increases,  while  that  of  Mg  is  nearly  constant  or  increases 
slightly,  with  decreasing  Cu  concentration.  Assigning  the 
effective  atomic  radius  of  element  i  as  half  the  mode  i-i  bond 
distance  (Table  III),  the  average  void  volume  and  packing 
fraction  near  i  atoms  were  calculated  for  each  alloy  as  the 
difference  between  the  respective  Voronoi  volume  and  the 
volume  of  atom  i,  and  as  the  fraction  of  the  atom  volume  rel¬ 
ative  the  VCV,  respectively  (Table  VI).  The  volume  of  voids 
and  the  packing  fraction  are  higher  near  larger  atoms.  The 
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FIG.  4.  Experimental  (dark  solid  lines)  and  QMD  simulated  (red  dashed 
lines)  RDFs  for  Ca60  xMg25Cui5+x  metallic  glasses.  The  curves  are  sepa 
rated  by  0.3  point  shifting  along  the  vertical  axis  for  better  view. 

void  volume  increases  but  the  packing  fraction  decreases 
near  Ca  and  Mg  atoms  with  a  decrease  in  Cu  concentration. 
The  packing  fraction  of  Cu  does  not  depend  on  the  alloy 
composition  within  the  25%-35%  Cu  range.  However,  it 
decreases  in  Ca60Mg25Cu15.  The  average  packing  fraction  of 
the  alloys  determined  from  their  density  (Table  VI)  shows 
similar  trends  and  increases  from  0.68  to  0.73  when  the  Cu 
content  increases  from  15%  to  35%.  These  values  show  that 
packing  is  significantly  more  efficient  than  in  monatomic, 
dense-random  packed  structures,  and  is  comparable  to  pack¬ 
ing  efficiencies  in  crystalline  metals. 

2.  Three-body  correlations 

Three-body  correlations  can  provide  additional  informa¬ 
tion  about  the  most  common  local  structures  around  a  cen¬ 
tered  atom  because  the  characteristic  angles  between  bonds 
connecting  any  two  atoms  located  in  the  first  or  second  coor¬ 
dination  shell  with  the  centered  atom  are  very  sensitive  to 


TABLE  III.  The  most  frequent  (mode)  distances  (in  Angstroms)  between 
pair  atoms  in  Ca  Mg  Cu  amorphous  alloys,  in  accord  to  QMD  and  CRMC 
simulations. 


Alloy 

Method 

Ca  Ca 

Ca  Mg 

Ca  Cu 

Mg  Mg 

Mg  Cu 

Cu  Cu 

Ca4oMg25Cu35 

QMD 

3.75 

3.42 

3.02 

3.14 

2.71 

2.47 

CRMC 

3.76 

3.45 

3.01 

3.14 

2.61 

2.49 

Ca45Mg25CU30 

QMD 

3.75 

3.41 

3.03 

3.13 

2.71 

2.48 

CRMC 

3.74 

3.43 

3.04 

3.12 

2.62 

2.49 

Ca50Mg25Cu25 

QMD 

3.75 

3.41 

3.07 

3.13 

2.75 

2.48 

CRMC 

3.75 

3.43 

3.04 

3.13 

2.73 

2.46 

Ca60Mg25Cu15 

QMD 

3.78 

3.46 

3.03 

3.08 

2.72 

2.51 

CRMC 

3.78 

3.47 

3.08 

3.12 

2.70 

2.50 

Average 

CRMC 

3.76 

3.45 

3.04 

3.12 

2.67 

2.49 

the  spatial  atom  distribution.22  Figure  9  illustrates  the  bond 
angle  distributions  in  the  Ca45Mg25Cu30  amorphous  alloy 
within  the  first  and  second  coordination  shells.  The  upper 
limits  of  the  bond  lengths  for  the  first  shell  were  set  to  the 
cut-off  values  given  in  Table  V,  while  the  bond  length  ranges 
for  the  second  shell  were  set  between  the  values  correspond¬ 
ing  to  the  first  and  second  minima  of  the  PRDFs  (see 
Figure  6).  Two  characteristic  maxima,  one  in  the  range  of 
50°  70°  and  another  between  100°  140°,  are  clearly  seen  in 
the  triplet  distributions  for  the  first  coordination  shell 
(Table  VII).  At  least  three  maxima  were  identified  in  bond 
angle  distributions  for  the  second  coordination  shell.  These 
are  near  30°,  60°,  and  90°  (Table  VIII). 

IV.  DISCUSSION 

The  atomic  structure  of  a  ternary  metallic  glass  is 
described  by  6  independent  PRDFs.  To  obtain  these  PRDFs 
analytically,  minimum  6  independent  scattering  experiments 
are  required.  However,  this  would  be  a  very  expensive  and 
time  consuming  research.  Moreover,  computer  modeling 
would  still  be  needed  to  produce  a  statistically  acceptable 


TABLE  II.  Partial  and  total  coordination  numbers  around  Ca,  Mg,  and  Cu  atoms  in  the  Ca  Mg  Cu  amorphous  alloys,  in  accord  to  QMD  and  RMC  simulations 
and  Voronoi  (Vor)  tessellation. 


CN 

Ca40Mg25Cu35 

Ca45Mg25Cu3o 

Ca5oMg25Cu25 

Ca60Mg25Cu15 

QMD 

CRMC 

Vor 

QMD 

CRMC 

Vor 

QMD 

CRMC 

Vor 

QMD 

CRMC 

Vor 

Ca  Ca 

7.0 

6.9 

7.0 

7.8 

7.1 

7.1 

8.6 

7.9 

7.9 

9.2 

8.6 

8.6 

Ca  Mg 

4.3 

3.7 

3.6 

4.1 

3.7 

3.8 

3.8 

3.4 

3.4 

3.5 

3.3 

3.4 

Ca  Cu 

4.9 

4.4 

4.4 

3.9 

3.8 

3.8 

3.2 

3.1 

3.1 

1.6 

1.7 

1.7 

Total 

16.2 

15.0 

15.0 

15.8 

14.6 

14.7 

15.6 

14.4 

14.4 

14.3 

13.6 

13.7 

Mg  Ca 

6.9 

5.9 

5.8 

7.3 

6.6 

6.9 

7.7 

6.8 

6.9 

8.5 

8.0 

8.0 

Mg  Mg 

2.6 

3.1 

3.0 

2.3 

2.8 

2.8 

2.1 

2.7 

2.8 

1.8 

2.7 

2.6 

Mg  Cu 

3.2 

3.7 

4.0 

2.5 

3.4 

3.3 

2.1 

3.0 

2.8 

1.1 

1.3 

1.6 

Total 

12.8 

12.7 

12.8 

12.2 

12.8 

13.0 

11.8 

12.5 

12.5 

11.3 

12.0 

12.2 

Cu  Ca 

5.6 

5.0 

5.1 

5.9 

5.6 

5.7 

6.4 

6.2 

6.2 

6.4 

6.8 

6.7 

Cu  Mg 

2.3 

2.7 

2.9 

2.1 

2.8 

2.7 

2.1 

3.0 

3.0 

1.8 

2.4 

2.6 

Cu  Cu 

1.6 

2.6 

2.1 

1.3 

1.9 

2.0 

1.0 

1.2 

1.5 

0.6 

0.6 

0.8 

Total 

9.4 

10.3 

10.1 

9.3 

10.4 

10.4 

9.5 

10.4 

10.7 

8.8 

9.8 

10.1 

6 

Distribution  Statement  A.  Approved  for  public  release;  distribution  unlimited. 


9wg-Cu(r)  3.  9Ca-Cu(0  SI  9ca-Ca(r)  3-  , 


1 2351 5-7  Senkov  et  at. 


J.  Appl.  Phys.  111,1 2351 5  (2012) 


FIG.  5.  Experimental  (dark  solid  lines)  and 
CRMC  simulated  (red  dashed  lines)  (a) 
total  structure  factors  and  (b)  RDFs  for 
Ca60  xMg25Cui5+x  metallic  glasses.  The 
curves  are  separated  by  shifting  along  the  verti 
cal  axis  for  better  view. 
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FIG.  6.  CRMC  simulated  (a)  Ca  Ca,  (b)  Ca 
Mg,  (c)  Ca  Cu,  (d)  Mg  Mg,  (e)  Mg  Cu,  and  (f) 
Cu  Cu  partial  radial  distribution  functions, 
gij(r)  for  Caso  xMg25Cu15+x  metallic  glasses. 
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FIG.  7.  Distribution  of  different  types  of  (a)  Ca  centered,  (b)  Mg  centered, 
and  (c)  Cu  centered  clusters  in  Ca6o  xMg2sCui5+x  amorphous  alloys. 


atomic  structure  of  the  glass.  Only  a  few  publications  on  bi¬ 
nary  systems  have  used  this  approach,  while  the  majority  of 
researchers  employ  proper  combinations  of  experiment  and 
simulation  to  determine  the  atomic  structure  of  metallic 
glasses  (for  review,  see  Ref.  22).  For  example,  it  has  been 
proven  in  a  number  of  publications  that  a  self-contained 
QMD  simulation,  which  requires  no  input  from  experiment 
and  thus  experiment-independent,  is  able  to  produce 
fairly  reasonable  (i.e.,  consistent  with  experimental  data) 
PRDFs.20'22  Thus  running  the  QMD  modeling  for  a  ternary 
alloy  to  obtain  6  independent  PRDFs  can  be  considered,  to 
some  extent,  to  be  equivalent  to  conducting  6  independent 
experiments.  Accordingly,  QMD  was  employed  in  our  work 
to  produce  6  independent  PRDFs,  consistent  with  the  experi¬ 
mental  data,  for  each  of  the  studied  Ca-Mg-Cu  amorphous 
alloys  (see  Figure  4).  Because  QMD  has  some  limitations, 
such  as  limited  number  of  atoms  and  high  quenching  rate, 
CRMC  was  further  used  to  produce  the  more  representative 
amorphous  structures,  which  agree  with  both  experimental 
and  QMD  results  The  use  of  a  synergistic  combination  of  the 
neutron  diffraction,  QMD  simulation,  and  CRMC  modeling 
allowed  us  to  resolve  the  amorphous  structures  of  several  ter¬ 
nary  Ca-Mg-Cu  alloys  in  the  most  efficient  manner,  thus  sav¬ 
ing  in  both  time  and  money,  which  is  in  line  with  a  currently 
widely  accepted  Integrated  Computational  Materials  Science 
and  Engineering  (ICMSE)  strategy.45  This  approach  pro¬ 
vided  consistent  results  for  different  alloy  compositions, 
which  additionally  supports/validates  the  appropriateness  of 
use  of  the  combination  of  experiment  and  modeling  in 
resolving  the  amorphous  structure.  It  is  necessary  to  point 
out  that  similar  approach  has  already  been  used  to  resolve 
the  amorphous  structure  of  binary  metallic  glasses.15'39'46  It 
is  used  for  the  first  time  in  this  paper  to  analyze  the  amor¬ 
phous  structure  of  ternary  metallic  glasses. 

QMD  simulated  structures  give  smaller  CNs  for  Cu-Cu, 
Cu-Mg,  Mg-Mg,  and  Mg-Cu  and  larger  Ca-Ca,  Ca-Mg, 
Ca-Cu,  Mg-Ca,  and  Cu-Ca  CNs  than  does  CRMC  (Table  II). 
QMD  simulated  structures  also  give  larger  mean  (Table  IV) 
and  cut-off  (Table  V)  pair  atom  distances.  Given  the  very 
high  QMD  quench  rate  and  the  periodic  boundary  condition 
on  a  relatively  small  volume,  these  minor  discrepancies  in 
structural  details  are  not  surprising.  At  the  same  time,  CNs 
obtained  with  Voronoi  tessellation  are  in  full  agreement  with 
CNs  obtained  directly  from  the  CRMC  partial  RDF’s.  Thus, 
only  the  CRMC-simulated  amorphous  structures  will  be  dis¬ 
cussed  here. 

A.  Pair  bond  distance  correlations 

rModc  values  (Table  III)  are  always  smaller  than  the  re¬ 
spective  rMean  values  (Table  IV).  The  difference  is  the  small¬ 
est  for  Ca-Ca  pairs  (3.2%)  and  increases  in  the  order  of  Ca- 
Mg  (4.3%),  Mg-Mg  (7.7%),  Ca-Cu  (8.2%),  Cu-Cu  (8.8%), 
and  Mg-Cu  (13.5%)  (Fig.  10(a)).  This  indicates  that  the  first- 
shell  gij(r)  peaks  are  asymmetric  (non-Gaussian),  especially 
for  Cu-containing  pairs.  rMode  and  rMea„  values  can  be  com¬ 
pared  with  atomic  separations  in  crystalline  metals  and 
alloys.  Crystalline  metallic  separations,  rM,  and  covalent  dis¬ 
tances,  rc,  were  estimated  as  a  sum  of  metallic  and  covalent 
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FIG.  8.  Fractions  of  3 , 4 , 5  ,  and  6  coordinated  vertices  in  Ca  ,  Mg  ,  and  Cu  centered  clusters  in  (a)  Ca40Mg25Cu35,  (b)  Ca45Mg25Cu3o,  (c)  Ca5oMg25Cu25, 
and  (d)  Ca6oMg25Cui5  amorphous  alloys. 


radii  of  respective  elements  (Table  IX).  The  metallic  radii 
are  half  the  nearest  distance  between  atoms  in  crystal  lattices 
of  the  respective  pure  metals.70  The  covalent  atomic  separa¬ 
tions71  were  obtained  from  analysis  of  nearest  distances  in 
crystalline  intermetallic  compounds  using  the  Cambridge 
Structural  Database.  From  Table  IX,  the  rMea„  and  val¬ 
ues  for  Ca-Ca  atom  pairs  are  smaller  than  rM,  so  that  the 
mean  and  mode  Ca-Ca  separations  in  Ca-Mg-Cu  glasses  are 
slightly  compressed.  However,  for  all  other  atomic  pairs, 
’'Mean  values  are  larger  and  rMode  values  are  smaller  than  rM 
(Fig.  10(b)).  The  larger  mean  separations  suggest  that  some 
Mg  and,  especially,  Cu  atoms  become  slightly  displaced 
from  each  other.  It  is  likely  that,  similar  to  Zn  in  Ca-Mg-Zn 
amorphous  alloys,23  some  Cu  atoms  occupy  gaps  between 
the  larger  Ca  and  Mg  atoms  in  the  first  coordination  shell 
and  their  distances  from  the  center  atom  vary  depending  on 
the  gap  sizes.  At  the  same  time,  rMode  f°r  the  atom  pairs  con¬ 
taining  Cu  are  even  smaller  than  the  respective  covalent 
bond  distances  in  crystalline  compounds,  and  the  mode  sepa¬ 
rations  of  Ca-Cu  and  Mg-Cu  are  well  represented  by  the 
comparison  between  rMode  and  rc  (Fig.  10(c)).  These  data 
show  that  the  absence  of  long-range  order  constraints  allows 
shortening  of  first  neighbor  interatomic  distances  in  amor¬ 
phous  structures,  as  compared  to  bond  distances  in  crystal¬ 
line  metals  and  alloys. 


The  Mg-Cu  (2.67  A)  and  Ca-Cu  (3.04  A)  mode  pair  dis¬ 
tances  are  much  shorter  than  those  calculated  from  the  Ca- 
Ca,  Mg-Mg,  and  Cu-Cu  rMode  values  (i.e.,  2.81  A  and  3.13  A, 
respectively,  see  Table  IX).  Shortening  of  Mg-Cu  and  Y-Cu 
bonds  has  also  been  observed  in  amorphous  Mg60Cu30Y  10  47 
Bond  shortening  has  been  discussed  in  amorphous  and  crys¬ 
talline  systems,48-50  and  is  likely  related  to  sp-d  electron 
hybridization.51  This  shortening  is  often  associated  with  a 
covalent  bonding  component  and  indicates  strong  chemical 
interactions,  suggesting  that  attractive  forces  will  increase  in 
the  sequence:  Ca-Mg,  Ca-Cu,  and  Mg-Cu.  It  has  been  sug¬ 
gested  that  the  thermodynamic  heat  of  mixing  can  represent 
this  chemical  interaction,52  and  experimental  data53-55  give 
the  heats  of  mixing  of  equimolar  binary  alloys 
as  —  6.0  ±  1 .0 kJ/mol  for  Ca-Mg,  —  6.5  ±  l.Okl/mol  for 
Ca-Cu,  and  —9.0  ±  1.0  kJ/mol  for  Cu-Mg.  These  values  give 
the  same  trend  as  for  bond  shortening,  but  the  link  between 
heats  of  mixing  and  bond  lengths  is  phenomenological.  Fur¬ 
ther,  the  heats  of  mixing  give  a  global  representation  of 
interatomic  bond  enthalpies,  which  include  the  number  of 
bonds  formed  and  the  energies  of  bonds  between  two  atoms. 
Both  of  these  values  are  likely  to  depend  upon  composition 
and  local  atom  environment,  so  that  trends  in  relative  bond 
lengths  obtained  from  heats  of  mixing  should  be  used  with 
caution. 
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Theta  (degrees)  Theta  (degrees) 

FIG.  9.  Ca  Ca  Ca,  Ca  Mg  Ca,  Ca  Cu  Ca,  Mg  Mg  Mg,  and  Cu  Cu  Cu  bond  angle  distributions  in  (a)  the  first  coordination  shell  and  (b)  second  coordination 
shell  in  the  Ca45Mg25Cu3o  amorphous  alloy. 


The  observation  that  nearest-neighbor  bond  lengths  are 
shorter  than  those  in  the  competing  long-range  ordered  crys¬ 
tals  may  help  explain  why  these  Ca-Mg-Cu  alloys  are  easy 
glass  formers.  The  equilibrium  crystalline  phases  that  com¬ 
pete  with  the  amorphous  phase  appear  to  minimize  total 
energy  by  sacrificing  optimal  short-range  order  to  achieve 
beneficial  long-range  order.  On  the  other  hand,  the  absence 
of  long-range  atomic  order  allows  metallic  glasses  to  reduce 
the  energy  difference  between  the  metastable  glass  and  equi¬ 
librium  crystal(s)  by  optimizing  short-range  atomic  interac¬ 
tions  and  arranging  the  atoms  in  efficiently  packed 
clusters.40,42,56'57  The  largest  contribution  to  condensed 
phase  stability  is  expected  to  come  from  nearest-neighbor 
interactions  and  volume  minimization,  so  that  the  optimized 
short-range  interactions,  as  well  as  the  presence  of  MRO,  in 
metallic  glasses  can  produce  structures  that  have  only  a  small 
energetic  disadvantage  relative  to  the  equilibrium  crystalline 
structure.22,23  Kinetic  constraints  from  quenching  restrict  the 
long  range  atomic  redistribution  needed  to  achieve  long- 
range  order  and  to  further  minimize  the  total  system  energy, 
thus  favoring  glass  formation.  The  present  results  suggest 
that  metallic  glasses  not  only  have  short  range  order  but  may 


TABLE  IV.  Weighted  average  (mean)  distances  (in  Angstroms)  between 
pair  atoms  in  the  first  shell  in  Ca  Mg  Cu  amorphous  alloys,  in  accord  to 
QMD  and  CRMC  simulations. 


Alloy 

Ca  Ca 

Ca  Mg 

Ca  Cu 

Mg  Mg 

Mg  Cu 

Cu  Cu 

Ca4oMg25Cu35 

QMD 

3.95 

3.64 

3.21 

3.31 

2.83 

2.59 

CRMC 

3.89 

3.55 

3.19 

3.34 

3.03 

2.85 

Ca45Mg25CU30 

QMD 

3.98 

3.63 

3.22 

3.28 

2.84 

2.59 

CRMC 

3.85 

3.61 

3.27 

3.37 

3.19 

2.75 

Ca5oMg25Cu25 

QMD 

4.01 

3.62 

3.25 

3.24 

2.89 

2.58 

CRMC 

3.89 

3.59 

3.33 

3.30 

3.15 

2.63 

Ca6oMg25Cu15 

QMD 

4.00 

3.64 

3.19 

3.29 

2.87 

2.60 

CRMC 

3.89 

3.66 

3.38 

3.46 

2.75 

2.60 

Average 

CRMC 

3.88 

3.60 

3.29 

3.37 

3.03 

2.71 

have  “better”  combination  of  short-range  topological  and 
chemical  order  (in  terms  of  optimal  bond  length  and/or 
atomic  arrangement)  than  the  competing  crystals. 

B.  Coordination  numbers  and  local  chemical  order 

Increasing  the  Cu  concentration  from  15  to  35  at.  %  con¬ 
tinuously  increases  the  total  coordination  number  around  Ca, 
CNca,  from  13.6  to  15.0  due  to  a  faster  increase  in  the  number 
of  smaller  Cu  atoms  {CNca-Cu  increases  from  1.7  to  4.4)  and  a 
slower  decrease  in  the  amount  of  the  larger  Ca  atoms  ( CNca-Ca 
decreases  from  8.6  to  6.9)  (Figure  1 1).  On  the  other  hand,  the 
CNMg  and  CNCu  almost  do  not  depend  on  alloy  composition. 
However,  similar  to  the  environment  around  Ca,  CN^g-cu  and 
CNcu-cu  increase  from  1.4  to  3.7  and  from  0.6  to  2.6,  respec¬ 
tively,  while  CNMg_Ca  and  CNcu-Ca  decrease  from  8.0  to  5.9 
and  from  6.8  to  5.0,  respectively,  with  the  increase  in  Cu 
concentration.  The  number  of  Mg  atoms  around  Ca,  Mg, 
and  Cu  very  weakly  depends  on  alloy  composition  and,  on  av¬ 
erage,  CNCa-Mg=  3.5  ±  0.2,  NMg.Mg  =  2.8  ±  0.3,  and  CNCu-Mg 
=  2.7  ±  0.3  (see  Table  II  and  Figure  1 1),  which  may  indicate  a 
rather  homogeneous  Mg  distribution.  The  dependences  of  the 
partial  coordination  numbers  on  Cu  concentration  can  be 


TABLE  V.  The  first  shell  cut  off  distances  between  pair  atoms  (in  Ang 
stroms)  in  Ca  Mg  Cu  amorphous  alloys,  in  accord  to  QMD  and  CRMC 
simulations. 


Alloy 

Ca  Ca 

Ca  Mg 

Ca  Cu 

Mg  Mg 

Mg  Cu 

Cu  Cu 

Ca4oMg25Cu35 

QMD 

4.88 

4.63 

4.26 

4.22 

3.58 

3.25 

CRMC 

4.74 

4.45 

4.11 

4.29 

4.04 

3.75 

Ca45Mg25Cu3o 

QMD 

4.91 

4.59 

4.26 

4.09 

3.58 

3.34 

CRMC 

4.73 

4.63 

4.31 

4.26 

4.34 

3.54 

Ca5oMg25CU25 

QMD 

4.99 

4.58 

4.40 

3.97 

3.84 

3.42 

CRMC 

4.82 

4.54 

4.44 

4.20 

4.30 

3.32 

Ca6oMg25Cu15 

QMD 

4.94 

4.58 

4.13 

4.05 

3.89 

3.31 

CRMC 

4.84 

4.70 

4.41 

4.46 

3.92 

3.22 
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TABLE  VI.  Voronoi  cell  volume,  void  volume,  and  packing  fraction  of  Ca,  Mg,  and  Cu  atoms  in  the  amorphous  structure  of  Ca  Mg  Cu  alloys.  The  packing 
fractions  of  the  alloys  are  also  given  in  the  last  column. 


Voronoi  volume  (A3) 

Void  volume  (A3) 

Packing  fraction 

Ca 

Mg 

Cu 

Ca 

Mg 

Cu 

Ca 

Mg 

Cu 

Alloy 

Ca40Mg25CU35 

35.99 

24.12 

13.50 

8.16 

6.64 

5.42 

0.78 

0.73 

0.61 

0.73 

Ca45Mg25Cll30 

37.19 

23.82 

13.86 

9.80 

7.92 

5.58 

0.74 

0.67 

0.60 

0.69 

Ca5oMg25Cu25 

38.71 

24.45 

13.16 

10.65 

9.30 

5.27 

0.73 

0.63 

0.61 

0.69 

Ca60Mg25Cui5 

40.06 

25.22 

14.47 

12.22 

9.32 

6.28 

0.70 

0.64 

0.57 

0.68 

described  by  linear  equations  (Table  X).  Extrapolation  to  low 
Cu-concentrations  predicts  no  Cu-Cu  first-neighbor  interac¬ 
tions  in  alloys  with  Cu  concentrations  less  than  10.3  at.  %.  At 
these  low  Cu  concentrations,  CNcu-Ca  and  CNcu-Mg  are  pre¬ 
dicted  to  be  7.3  and  2.5,  respectively,  CNra  =  13  and 
CNMg  =  11.8. 

To  determine  the  degree  of  chemical  short  range  order 
(CSRO)  around  Ca,  Mg,  and  Cu  atoms,  the  Warren-Cowley 
parameter,  ay,  which  is  defined  as58 

a,7  =  1  CNy/ (cjCNj)  (7) 

was  used.  Here,  cj  is  the  atomic  fraction  of  the  element  j, 
CNij  is  the  partial  coordination  number  of  element  j  around 
the  element  i,  and  CNt  is  the  total  coordination  number 
around  element  i.  A  negative  ay  (/  #  j)  indicates  the  pres¬ 
ence  of  CSRO  (i.e.,  the  number  of  atoms  j  in  the  first  shell 
exceeds  the  average  concentration).  For  a  random  solution, 
a jj  are  zero.  A  modified  CSRO  parameter,  a,^,  which  is 
defined  as 

«,■(/*)  =  1  {CNij  +  CNik)/  ({cj  +  ct)OV/)  for  j  ±  k 

(8) 

can  also  be  used  to  analyze  the  CSRO  in  ternary  alloys.  Neg¬ 
ative  rJ.j(jk)  values  indicate  the  presence  of  CSRO,  while  posi¬ 
tive  values  indicate  the  presence  of  chemical  short  range 
clustering  (CSRC)  near  the  i  atom.  7.!(jkl  and  a„  are  dependent 
parameters,  namely 

(1  -Ci)ai(jk)  =  -Clan.  (9) 

Therefore,  a„  <  0  should  indicate  CSRC  (i.e.,  increased  local 
concentration  of  like  atoms)  and  a„  >  0  should  indicate 
CSRO  around  atom  i. 

Values  of  ay  and  ai(jk)  for  Ca-Mg-Cu  glasses  are  given 
in  Tables  XI  and  XII,  respectively.  The  alloys  show  CSRC 
of  Ca  atoms  that  is  largest  for  the  alloy  with  40%  Ca 
(aCaca=  —0.16)  and  it  tends  to  decrease  with  an  increase  in 
Ca  concentration  (at  60%  Ca  aCaCa  =  —  0.05).  Pronounced 


CSRO  is  detected  near  Cu  atoms.  These  local  regions  are 
enriched  with  Ca  and  Mg,  supporting  the  earlier  observation 
of  strong  Ca-Cu  and  Mg-Cu  bonding.  Near  Mg  atoms,  aMgCa 
is  negative  and  aMgCu  is  positive,  indicating  enrichment  in 
Ca  and  a  Cu  deficit  in  Mg-centered  clusters.  The  ai(jk)  crite¬ 
rion  shows  clear  evidence  of  CSRC  near  Ca  atoms  and 
CSRO  near  Cu  atoms  when  unlike  elements  are  considered 
together  (Table  XII).  At  the  same  time,  ai(jk)  is  almost  zero 
or  slightly  negative  for  Mg-centered  clusters  indicating  neu¬ 
tral  environment. 

C.  Voronoi  tessellation  analysis 

Voronoi  analysis  indicates  that  Kasper-type  polyhedra 
dominate  in  Ca-Mg-Cu  glass  structures.  Among  the  clusters 
in  Figure  7,  (0, 4,4,0),  (0,3, 6,0),  (0, 2,8,0),  (0,2, 8,1), 
(0,0,12,0),  and  (0,1,10,2)  are  non-distorted  Kasper  polyhe¬ 
dra,  whereas  (0,3,6, 1),  (0, 4,4,2),  (0,3, 6, 2),  (0,4, 4,3), 
(0,2, 8, 2),  (0,3, 6, 3),  (0,4,4, 4),  (0,2, 8,3),  (0,3, 6, 4),  (0,4, 4,5), 
(0,1,10,3),  (0,2, 8, 4),  (0, 3,6,5),  (0, 4,4,6),  (0,1,10,4),  (0,2, 8,5), 
and  (0,3, 6, 6)  are  distorted  Kasper  polyhedra  with  four-fold 
and  six-fold  disclinations,22  and  only  (0,3, 7, 4),  (1,2, 6, 3), 
(1,2, 5,2),  (1,3, 3, 2),  and  (1,2, 5, 3)  are  other-type  polyhedra. 
The  Kasper  polyhedra  edges  can  have  different  lengths  this 
is  inevitable  in  glass;  therefore,  they  are  not  identical  even 
though  they  may  share  the  same  topology  (in  terms  of  the 
Voronoi  index).  Similar  to  other  work,22  a  cluster  is  identi¬ 
fied  here  as  a  non-distorted  Kasper  polyhedron  if  it  has  a 
Voronoi  index  of  the  Kasper  polyhedron  and  the  term 
“distorted”  is  assigned  to  Kasper  polyhedra  with  extrinsic 
disclinations. 

Kasper  polyhedra  and  their  distorted  variants  account 
for  over  50%  of  the  nearest-neighbor  clusters  (Table  XIII). 
Since  Kasper  polyhedra  are  polytetrahedral,  and  many  other 
polyhedra  also  contain  tetrahedra,  we  conclude  that  the  topo¬ 
logical  SRO  of  Ca-Mg-Cu  is  polytetrahedral  in  nature,  which 
has  been  demonstrated  in  many  other  metallic  glasses.22  The 
current  study,  therefore,  further  supports  the  idea  that  the 
SRO  of  MGs  is  characterized  by  polytetrahedral  packing  via 
Kasper  clusters  and  their  distorted  variants.  The  regularity  of 


TABLE  VII.  Positions  of  the  first  and  the  second  maxima  (in  degrees)  in  the  bond  angle  distributions  of  five  triplets  within  the  first  coordination  shell  in  four 
Ca  Mg  Cu  amorphous  alloys. 


1st  shell 


Ca  Ca  Ca 


Ca  Mg  Ca 


Ca  Cu  Ca 


Mg  Mg  Mg 


Cu  Cu  Cu 


Ca40Mg25CU35 

57.0,  101.1 

62.5,  111.1 

71.7,  122.6 

57.4,  108.2 

52.8,  110.5 

Ca45Mg25Cu3o 

58.0,  103.8 

62.3,  111.1 

69.6,  121.4 

57.4,  110.4 

55.4,  111.0 

Ca5oMg25Cu25 

58.2,  104.3 

62.5,  113.1 

69.8,  122.3 

57.2,  108.5 

55.5,  110.5 

Ca60Mg25Cu15 

58.4,  105.5 

62.3,  113.3 

66.5,  121.6 

55.1.  108.9 

61.2,  107.9 
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TABLE  VIII.  Positions  of  three  maxima  (in  degrees)  observed  in  the  bond  angle  distributions  of  five  triplets  within  the  second  coordination  shell  in  four  Ca 
Mg  Cu  amorphous  alloys. 


2nd  shell 

Ca  Ca  Ca 

Ca  Mg  Ca 

Ca  Cu  Ca 

Mg  Mg  Mg 

Cu  Cu  Cu 

Ca4oMg25Cu35 

Ca45Mg25Cu3o 

Ca5oMg25Cu25 

Ca60Mg25Cui5 

29.5,  62.5,  88.5 
30.5,61.0,  85.0 

30.5,  60.0,  86.5 

30.5,  60.5,  84.5 

31.5,63.0,88.5 
31.5,60.0,  92.5 
31.5,60.0, 93.5 
31.5,61.5,90.0 

33.5.61.5,  86.5 

33.5,  60.0,  90.0 

32.5.62.5,  87.5 
31.5,62.5,94.5 

27.5,  63.5,  86.5 

27.5,  60.0,  90.0 

27.5,  60.0,  90.5 

26.5,  57.5,  85.0 

24.5,  60.0,  90.0 

25.5,  57.5,  90.0 

23.5,  60.0,  84.5 
24.5,58.5,95.5 

polytetrahedral  packing  can  be  evaluated  using  the  ratio 
between  the  number  of  regular  Kasper  polyhedra  and  total 
number  of  regular  and  distorted  Kasper  polyhedra.  We  find 
that  the  atomic  packing  around  smaller  atoms  (Cu)  is  much 
more  regular  polytetrahedral  than  around  larger  atoms  (Ca). 
For  example,  the  fractions  of  regular  Kasper  polyhedra 
around  Cu  is  about  55%,  around  Mg  range  from  23%  to 
32%,  and  decreases  to  13%-20%  for  Ca-centered  clusters. 
This  is  reminiscent  of  Cu-Zr,  where  the  packing  around  Cu 
appears  much  more  regular  than  the  packing  around  Zr.22 
The  more  regular  packing  around  Cu  can  indicate  that 
Cu-centered  clusters  are  the  primary  structure-forming  clus¬ 
ters,  whereas  Mg-  and  Ca-centered  clusters  are  secondary 
clusters  resulting  from  specific  arrangements  of  atoms  in  the 
vertices  of  and  voids  between  the  Cu  clusters  (see  Sec.  IV  E 
below).  This  can  also  be  explained  by  very  strong  interac¬ 
tions  of  Cu  with  Ca  and  Mg  resulting  in  bond  shortening. 
This  feature  could  be  quite  general  for  metallic  glasses,  i.e., 
polytetrahedral  packing  is  more  prominent  and  better  estab¬ 
lished  around  smaller  atoms.  This  is  probably  because 
smaller  atoms  have  smaller  CNs,  and  thus  have  limited 
choices  for  efficient  packing,  in  comparison  to  larger  atoms 
with  more  neighbors  and  possibilities. 

Similar  to  Ca-Mg-Zn  BMGs,2  the  amorphous  structures 
of  Ca-Mg-Cu  glasses  contain  a  very  low  fraction  of  icosahe- 
dral  SRO,  in  spite  of  the  observation  that  five-coordinated  ver¬ 
tices  dominate  in  all  clusters.  A  large  fraction  of  five-fold 
bonds  and  the  lack  of  icosahedral  SRO  have  also  been  found 
in  the  MD-simulated  amorphous  structure  of  Mg-Cu  alloys.24 
The  dominance  of  five-fold  vertices  suggests  that  tetrahedra 
prefer  to  cluster  into  pentagonal  bi-pyramids.  Analysis  of 
amorphous  Mg60Cu30Y  10  (Ref.  47)  shows  a  prevalence  of 
fivefold  bonds,  with  local  atom  arrangements  very  similar  to 
the  competing  Mg2Cu  (Laves  phase)  and  YCu2  crystals.  Simi- 

TABLE  IX.  Metallic,  rm,10  and  covalent,  /‘c,71  crystalline  bond  distances 
between  Ca,  Mg,  and  Cu  atom  pairs.  These  distances  are  compared  to  rMode 
and  /'Mean  from  the  CRMC  PRDFs.  rMode  and  rMearl  are  the  average  values 
for  the  four  studied  alloys. 


Ca  Ca 

Ca  Mg 

Ca  Cu 

Mg  Mg 

Mg  Cu 

Cu  Cu 

rM  (A) 

3.94 

3.57 

3.25 

3.20 

2.88 

2.56 

rc  (A) 

3.52 

3.17 

3.08 

2.82 

2.73 

2.64 

^Mode  (A) 

3.76 

3.45 

3.04 

3.13 

2.67 

2.49 

^Mean  (A) 

3.88 

3.60 

3.29 

3.37 

3.03 

2.71 

100%  (rMean/^Mode  1) 

3.2 

4.3 

8.2 

7.7 

13.5 

8.8 

100%  (rM/rMode  1) 

4.8 

3.5 

6.9 

2.2 

7.9 

2.8 

100%  (rc/rMode  1) 

6.4 

8.1 

1.3 

9.9 

2.2 

6.0 

100%  (rM/rMean  1) 

1.5 

0.8 

1.2 

5.0 

5.0 

5.5 

lar  comparison  of  coordination  polyhedra  in  the  amorphous 
structure  of  Ca60MgxZn4o  x  with  polyhedra  in  the  competing 
crystal  phases,  CaMg2  and  CaZn,59  has  shown  that  SRO 
arrangements  around  Ca  and  Mg  atoms  were  different  from 
those  in  the  crystal  phases.23  However,  the  (0, 3,6,0)  Zn- 
centered  coordination  polyhedron,  typical  to  CaZn,  was  one 
of  the  most  common  in  the  amorphous  structure. 

Amorphization  of  Ca6oMg25Cu15  competes  with  CaMg2 
compound  formation,  amorphization  of  Ca40Mg25Cu35 
competes  with  the  Cu2Mg  phase,  while  solidification  of 
Ca45Mg25Cu3o  and  Ca5oMg25Cu25  alloys  forms  three  equally 
probable  crystal  phases,  CaMg2  Cu2Mg  and  CaCu.5’6  CaMg2 
is  a  hexagonal  Laves  phase  (space  group  P63/mmc,  space 
group  number  194)  with  lattice  parameters  a  =5. 170  A  and 
c  =  8.50  A.  This  crystal  structure  contains  three  coordination 
polyhedra,  two  of  which  are  Mg-centered  (0,0,12,0)  icosahe- 
dra  with  6  Mg  and  6  Ca  atoms  in  the  first  shell,  and  the  third 
is  a  Ca-centered  (0,0,12,4)  Friauf  polyhedron  containing 
12  Mg  and  4  Ca  atoms  in  the  first  shell.  Cu2Mg  is  a  cubic 
Laves  phase  (space  group  is  Fd-3  m,  space  group  number 
227)  with  the  lattice  parameter  a  =  6.990  A.  The  crystal  struc¬ 
ture  of  this  phase  has  two  coordination  polyhedra:  a  Cu- 
centered  (0,0,12,0)  icosahedron  with  6  Cu  and  6  Mg  atoms  in 
the  first  shell,  and  a  Mg-centered  (0,0,12,4)  Friauf  polyhedron 
with  12  Cu  and  4  Mg  atoms  in  the  first  shell.  The  CaCu  phase 
has  a  primitive  crystal  structure  (space  group  P2!/m,  space 
group  number  11)  with  lattice  parameters  a  =19.47  A, 
b  =  4.271  A,  c  =  5.880  A,  and  /?  =  94.30°.  It  has  10  character¬ 
istic  coordination  polyhedra,  5  of  which  are  Cu-centered  Kas¬ 
per  polyhedra  (0,3, 6,0)  with  2  Cu  and  7  Ca  atoms  in  the  first 
shell,  and  the  other  5  are  Ca-centered  polyhedra  (3, 1,6, 6,1), 
(6, 1,2,2, 4, 0,2),  (1,3, 6, 4, 3),  (3,3, 3, 4, 3,0,1),  and  (1,3, 7,6),  with 
7  Cu  and  10  Ca  atoms  in  the  first  shell.60  (Additional  indices 
in  the  Voronoi  signatures  of  the  Ca-centered  polyhedra  given 
above  correspond  to  u7,  v$,  and  v9,  respectively.) 

None  of  the  Ca-centered  coordination  polyhedra  of  the 
competing  crystal  phases  are  present  in  amorphous  structures 
of  Ca-Mg-Cu.  The  icosahedral  SRO  in  the  CaMg2  and  Cu2Mg 
crystal  phases  near  Mg  and  Cu  atoms,  respectively,  is  not  typi¬ 
cal  of  Ca-Mg-Cu  glass  structures.  The  (0,3, 6,0)  Cu-centered 
coordination  polyhedron  in  the  CaCu  phase  is  also  common  in 
the  studied  amorphous  structures.  This  (0,3 ,6,0)  Cu-centered 
cluster  is  probably  the  only  link  between  the  competing  crys¬ 
tal  and  amorphous  states  and,  therefore,  a  reduced  fraction  of 
this  cluster  in  the  amorphous  structures  should  indicate  better 
GFA.  Indeed,  the  Ca60Mg25Cu15  and  Ca4oMg25Cu35  amor¬ 
phous  alloys  have  much  higher  fraction  of  the  (0,3 ,6,0)  Cu- 
centered  clusters  than  the  better  glass  forming  Ca45Mg25Cu3o 
and  Ca5oMg25Cu25  alloys. 
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FIG.  10.  Comparison  of  the  mode  atomic  pair  distances  with  (a)  the  respective  mean  distances  in  the  amorphous  structure  of  the  studied  Ca  Mg  Cu  alloys,  and 
with  the  reported  (b)  metallic,  rM,  and  (c)  covalent,  rc,  bond  distances. 


D.  Three-body  correlations 

The  Voronoi  analysis  above  indicates  that  tetrahedra 
are  dominant  building  blocks  in  amorphous  structures  of 
Ca-Mg-Cu  alloys.  With  more  than  50%  frequency,  the  tetra¬ 
hedra  form  pentagonal  bi-pyramids  and  Kasper  polyhedra. 
Because  the  characteristic  vertex  angles  for  regular  tetrahe¬ 
dra  are  60°  and  for  the  regular  pentagonal  bi-pyramids 
are  60°  and  108°,  the  deviations  of  angle  peak  positions  in 
three-body  correlations  indicate  distortions  from  the  regular 
configurations.  These  distortions  are  mainly  caused  by  dif¬ 
ferent  radii  of  Ca,  Mg,  and  Cu  atoms,  as  well  as  non- 
Gaussian  distributions  of  the  bond  distances  in  the  first  coor¬ 
dination  shell. 

Table  VII  gives  the  first  and  second  maxima  in  bond 
angle  distributions  of  Ca-Ca-Ca,  Ca-Mg-Ca,  Ca-Cu-Ca, 
Mg-Mg-Mg,  and  Cu-Cu-Cu  triplets  in  the  first  coordination 
shell.  If  like  atoms  are  packed  closely,  then  the  characteristic 
angle  for  the  i-i-i  triplets  will  be  equal  or  close  to  60°.  The 
first  peak  for  Ca-Ca-Ca,  Mg-Mg-Mg,  and  Cu-Cu-Cu  triplets 


occurs  at  —57.9°,  56.8°,  and  56.2°,  respectively.  The  peak 
shift  to  angles  less  than  60°  indicates  that  the  i-i  distances 
within  the  first  shell  are  longer  than  the  bonds  between  the 
center  and  two  first-shell  i  atoms.  From  Sec.  IV  A,  we  take 
the  distances  between  the  center  and  shell  like  atoms  as  rMode 
and  the  distance  between  the  shell  atom  pair  as  rMean,  so  that 
the  angle  of  an  i-i-i  triplet  is  estimated  as 

(X  —  COS  (/  Mean /2t  Mode)  •  (10) 

The  calculated  angles  for  Ca-,  Mg-,  and  Cu-triplets  using 
averaged  rMode  and  rMean  values  from  CRMC  simulations 
(see  Table  IX)  are  58.9°,  57.4°,  and  57.0°,  respectively.  These 
are  all  within  1°  of  measured  angles,  giving  good  agreement. 

For  Ca-Mg-Ca  and  Ca-Cu-Ca  triplets,  the  first  peaks  are 
at  —62.4°  and  69.4°,  respectively.  These  larger  angles  result 
from  the  smaller  atomic  radii  of  Mg  and  Cu  atoms.  Using 
r Mode  f°r  Ca-Ca,  Ca-Mg,  and  Ca-Cu  separations,  the  Ca-Mg- 
Ca  and  Ca-Cu-Ca  triplet  angles  are  estimated  as  63.0°  and 
67.5°,  respectively.  This  is  very  good  agreement  for  Ca-Mg- 


FIG.  11.  Dependence  on  the  Cu  con 
centration  of  (a)  the  total  coordination 
number,  CNm,  where  M  is  Ca,  Mg,  or 
Cu  center  atom,  and  partial  coordina 
tion  numbers  (b)  CNM  Ca,  (c)  CNM  Mg, 
and  (d)  CNM  cu- 
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TABLE  X.  Linear  fits  to  the  dependence  on  Cu  content  (at.  %)  of  coordina 
tion  number  data  (Figure  11).  R2  is  the  coefficient  of  determination  of  the  re 
spective  linear  regression. 


Atomic  pair 

Gradient 

Intercept 

R2 

Ca  Ca 

0.090 

9.98 

0.963 

Ca  Mg 

0.022 

2.95 

0.884 

Ca  Cu 

0.125 

0.00 

0.992 

Mg  Ca 

0.101 

9.48 

0.978 

Mg  Mg 

0.018 

2.36 

0.639 

Mg  Cu 

0.110 

0.00 

0.956 

Cu  Ca 

0.089 

8.24 

0.966 

Cu  Mg 

0.0154 

2.32 

0.278 

Cu  Cu 

0.099 

1.02 

0.951 

Ca  and  shows  reasonable  agreement  for  Ca-Cu-Ca  given  the 
uncertainty  in  bond  angles. 

The  second  bond-angle  distribution  peak  for  Ca-Ca-Ca 
triplets  is  located  near  101°  in  the  Ca4oMg25Cu3g  alloy  and 
shifts  to  ~  1 06°  in  Ca60Mg25Cu15.  This  is  a  little  smaller  than 
the  108°  interior  angle  of  a  regular  pentagon  and  may  be  due 
to  distortion  when  some  vertices  are  occupied  by  smaller  Mg 
and  Cu  atoms.  This  also  explains  the  angle  decrease 
with  increasing  Cu  concentration.  Due  to  shorter  Ca-Cu  and 
Ca-Mg  bonds,  the  interior  angle  at  the  Cu  or  Mg  vertex 
increases  at  the  expense  of  other  interior  angles,  the  sum  of 
which  is  constant  at  540°.  Indeed,  the  second  peak  for  Ca- 
Mg-Ca  and  Ca-Cu-Ca  triplets  is  located  at  1 1 1°  —  1 14°  and  at 
121°  — 130°,  respectively.  Using  these  values  of  interior 
angles  at  Ca,  Mg,  and  Cu  vertices,  on  average  2.6  Ca, 
1.1  Mg,  and  1.3  Cu  atoms  in  the  Ca40Mg25Cu35  alloy  and 
~4.0  Ca,  0.7  Mg,  and  0.3  Cu  atoms  in  the  Ca60Mg25CuI5 
alloy  are  estimated  per  pentagon  to  satisfy  the  sum  of  vertex 
angles  of  540°. 

The  first  three  maxima  in  bond  angle  distributions  of 
Ca-Ca-Ca,  Ca-Mg-Ca,  Ca-Cu-Ca,  Mg-Mg-Mg,  and  Cu-Cu-Cu 
triplets  in  the  second  coordination  shell  are  given  in  Table  VIII 
and  Figure  9(b).  If  like  atoms  contact  in  the  second  shell,  then 
the  characteristic  first  angle  for  i-i-i  triplets  will  be  equal  or 
close  to  30°.  This  angle  is  estimated  for  Ca-,  Mg-,  and  Cu  trip¬ 
lets  by  assuming  the  i-i  pair  in  the  2nd  coordination  shell  are 
separated  by  rmode  and  that  the  second  shell  is  displaced  from 
the  central  i  atom  by  a  distance  ( I'Mode  +  ?),  where  r  is  the 
weighted  average  of  rmode  values  for  the  atoms  in  the  first  shell 
of  the  i  atom.  The  triplet  angle  is  estimated  as 

p  =  COS  '{1  -  ( 1  /2)  [r Mode  /  (r Mode  +  Clf}-  (H) 

The  angles  for  a  given  i-i-i  triplet  are  averaged  over  all  four 
compositions.  The  estimated  Ca-,  Mg-,  and  Cu  triplet  angles 
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TABLE  XII.  The  short  range  order  parameters,  ai(jk),  for  Ca  ,  Mg  ,  and 
Cu  centered  clusters  in  the  Ca  Mg  Cu  amorphous  alloys. 


ai(jk>  * 

Ca 

Mg 

Cu 

Ca40Mg25CU35 

0.11 

0.00 

0.15 

Ca45Mg25Cu3o 

0.06 

0.05 

0.16 

Ca50Mg25Cu25 

0.10 

0.05 

0.14 

Ca60Mg25Cui5 

0.08 

0.04 

0.10 

are  30.8°,  27.9°,  and  24.4°,  respectively,  and  are  all  within 
1°  of  the  measured  angles  (30.3°,  27.3°,  and  24.5°).  A  simi¬ 
lar  topological  analysis  of  Ca-Mg-Ca  and  Ca-Cu-Ca  triplets 
for  the  2nd  coordination  shell  give  estimated  angles  of  32.2° 
and  34.0°,  which  compare  well  with  the  measured  values  of 
31.5°  and  32.8°. 

A  second  peak  is  expected  near  60°  when  a  full  atom 
occurs  between  the  i-i  pair  in  the  second  shell,  and  the  bond 
angle  distributions  show  a  rather  broad  peak  centered  very 
near  60°  for  the  five  triplets  studied  (Figure  9(b)).  To  model 
this,  an  atom  with  the  effective  radius  of  atoms  in  the  1st 
shell  of  i  atoms  is  placed  between  the  i-i  atom  pair  in  the  2nd 
shell  so  that, 

P  =  cos  X{1  -  (l/2)[(rUode  +  r)/(rMode  +  f)]2}-  (12) 

By  inspection,  this  gives  an  angle  of  60°  for  each  i-i-i  triplet, 
in  agreement  with  angles  of  61.0°,  60.3°,  and  59.0°  from 
CRMC  for  Ca-,  Mg-,  and  Cu  triplets.  The  Ca-Mg-Ca  (61.1°) 
and  Ca-Cu-Ca  (61.6°)  angles  are  both  slightly  larger  than 
60°.  While  this  trend  is  matched  using  Eq.  (12)  for  these  trip¬ 
lets,  the  estimated  angles  are  somewhat  larger  (63.1°  and 
67.6°,  respectively). 

A  third  peak  in  the  second  coordination  shell  may  be 
expected  near  90°  by  summing  the  1st  and  2nd  peaks.  While  a 
local  maxima  are  shown  in  all  five  bond  angle  distributions,  it 
is  a  very  shallow  and  broad  distribution  of  angles.  The  many 
configurations  for  intervening  atoms  make  it  difficult  to  per¬ 
form  a  simple  topological  analysis  of  the  included  angles. 

E.  Medium  range  order 

The  neutron  diffraction  pre-peak  at  Q  «  1 .2  A  1  of  Ca- 
Mg-Cu  glasses  (Figure  1)  corresponds  to  certain  medium- 
range  correlations  in  real  space.22  Generally,  MRO  is  seen  in 
solute-lean  glasses15’61-64  and  is  explained  in  terms  of  a 
sublattice-like  pattern  formed  by  solutes42  and  significantly 
higher  scattering  amplitude  for  solute  than  for  solvent 
atoms.22  However,  similar  to  this  work,  a  pre-peak  has  also 
been  seen  in  concentrated  glasses,  such  as  CeggAl.jg  (Ref.  65) 
and  Ca60MgxZn40  x.2°  Our  analysis,  given  below  for  the 


TABLE  XI.  Short  range  order  parameters,  a^,  for  Ca  ,  Mg  ,  and  Cu  centered  clusters  in  the  Ca  Mg  Cu  amorphous  alloys. 


“U 

Ca  Ca 

Ca  Mg 

Ca  Cu 

Mg  Ca 

Mg  Mg 

Mg  Cu 

Cu  Ca 

Cu  Mg 

Cu  Cu 

Ca40Mg25Cu35 

0.16 

0.03 

0.16 

0.15 

0.00 

0.17 

0.21 

0.05 

0.28 

Ca45Mg25Cu3o 

0.07 

0.04 

0.15 

0.17 

0.15 

0.14 

0.21 

0.08 

0.38 

Ca5oMg25Cu25 

0.10 

0.05 

0.14 

0.09 

0.14 

0.05 

0.15 

0.12 

0.42 

Ca6oMg2sCu15 

0.05 

0.02 

0.17 

0.11 

0.11 

0.25 

0.17 

0.07 

0.56 
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TABLE  XIII.  Fractions  of  Ca  ,  Mg  ,  and  Cu  centered  Kasper  polyhedra,  distorted  Kasper  polyhedra,  and  sum  of  both. 


Alloy 

Kasper  polyhedra 

Distorted  Kasper  polyhedra 

Sum 

Ca 

Mg 

Cu 

Ca 

Mg 

Cu 

Ca 

Mg 

Cu 

Ca4oMg25CU35 

0.059 

0.142 

0.323 

0.408 

0.467 

0.275 

0.467 

0.610 

0.598 

Ca45Mg25CU30 

0.079 

0.152 

0.323 

0.424 

0.402 

0.269 

0.503 

0.555 

0.592 

Ca5oMg25Cu25 

0.102 

0.192 

0.308 

0.454 

0.425 

0.260 

0.556 

0.617 

0.568 

Ca6oMg25Cu15 

0.121 

0.190 

0.338 

0.498 

0.407 

0.283 

0.619 

0.597 

0.621 

Ca45Mg25Cu3o  metallic  glass,  shows  that  the  MRO  in  Ca- 
Mg-Cu  metal  glasses  can  be  explained  as  produced  by  a  face 
center  cubic  (FCC)  like  local  ordering  of  Cu-centered  clus¬ 
ters,  and  the  pre-peak  results  from  this  FCC-like  ordering 
combined  with  the  strong  neutron  scattering  from  Cu. 

Figure  12  shows  six  CRMC-simulated  partial  structure 
factors,  A,j(Q),  for  the  Ca45Mg25Cu3o  metallic  glass.  The 
dashed  vertical  line  in  these  figures  corresponds  to  the  posi¬ 
tion  of  the  pre-peak  maximum  (Q  =  1.25  A  1 )  on  the  respec¬ 
tive  total  S(Q)  (see  Figure  1(a)).  It  can  be  clearly  seen  that 
this  pre-peak  is  originated  from  strong  Cu-Cu  and  Mg-Cu 
correlations  (Figures  12(e)  and  12(f)),  as  well  as  a  weak  Mg- 
Mg  correlation  (Figure  12(d)).  No  peak  intensity  is  seen  in 
this  Q  range  for  Ca-Ca,  Ca-Mg,  and  Ca-Cu  correlations.  This 
result  indicates  that,  similar  to  solute-lean  metallic  glasses, 
MRO  in  Ca-Mg-Cu  metal  glasses  is  caused  by  solute-solute 
(mainly  Cu-Cu  and  Mg-Cu)  interactions. 

Following  a  recent  model,43  we  assume  that  MRO  is 
caused  by  local  arrangement  of  Cu-centered  clusters  in  an 
FCC  motif,  as  shown  schematically  in  Figure  13.  This  is 
supported  not  only  by  diffraction  results  but  also  by  the  tri¬ 
plet  angle  distributions  for  the  2nd  coordination  shell, 


which  show  a  surprising  preference  for  angles  of  60°  and 
90°  irrespective  of  composition.  An  FCC  cell  consists  of  4 
Cu-centered  clusters  with  4  octahedral  (/?)  and  8  tetrahedral 
(y)  interstitial  sites,  which  tentatively  can  also  be  occupied 
by  alloying  elements.  On  average,  each  Cu-centered  cluster 
consists  of  11.4  atoms,  so  that  the  number  of  atoms  per  unit 
cell  is  estimated  to  be  between  45.6  (no  occupied  interstitial 
sites)  and  57.6  (all  12  sites  are  occupied).  To  satisfy  the 
alloy  density  (pa=  0.0373  A  3)  the  cell  volume,  Vc,  should, 
therefore,  be  between  1223  A3  and  1544  A3,  or  the  cluster 
unit  cell  parameter,  ac,  should  be  between  10.7  A  and 
11.6  A.  This  results  in  the  closest  average  distance  between 
the  Cu-centered  clusters,  rCu-cuCluster  (=  aJy/2)  to  be 
between  7.6  A  and  8.2  A.  This  distance  corresponds  to  the 
third  peak  on  the  Cu-Cu  PRDF,  located  at  r  =  7.8  A  (see 
Figure  6(f)).  Assigning  this  value  to  rCu_cucluster,  the  param¬ 
eter  ac  is  refined  as  ac=  11.0  A,  which,  at  the  given  alloy 
density,  gives  49.6  atoms  per  super-cell.  This  means  that 
only  4  of  the  12  interstitial  sites  (i.e.,  ~33%)  per  unit  cell 
are  occupied.  Conducting  similar  analysis  for  the  Mg- 
centered  clusters  gives  the  average  distance  between  these 
clusters,  rMg_Mgcluster,  in  the  range  of  8.1  A  to  8.6  A.  No 


FIG.  12.  Partial  structure  factors,  Ajj(Q),  for  Ca45Mg25Cu3o  metallic  glass. 
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FIG.  13.  Schematic  of  a  face  centered  cubic  cell  formed  by  Cu  centered 
clusters  (shown  by  red).  Positions  of  octahedral  (/?)  and  tetrahedral  (y)  inter 
stitial  voids  are  also  shown. 

peak  is  present  on  the  Mg-Mg  PRDF  in  this  r  range  (see 
Figure  6(d)).  We  conclude  that  Mg-Mg  clusters  are  not 
arranged  in  an  FCC  motif. 

While  the  third  peak  on  the  Cu-Cu  PRDF  can  be 
assigned  to  rCu_Cucluster,  which  is  the  distance  between  the 
Cu-centered  clusters  in  the  j  0,0,0}  and  {1/2, 1/2,0)  positions 
of  the  FCC  super-cell,  the  wide  second  peak  on  the  Cu-Cu 
PRDF  can  be  assigned  to  the  distances  between  the  Cu  in 
the  center  of  the  cluster  and  Cu  located  in  the  nearest  tetra¬ 
hedral  and/or  octahedral  sites,  as  well  as  to  the  distance 
between  nearest  tetrahedral  and  octahedral  sites.  The  small 
fourth  Cu-Cu  PRDF  peak  located  at  r  =  9.2  A  can  be 
assigned  to  the  distance  between  a  tetrahedral  site  and  sec¬ 
ond  nearest  octahedral  site  or  second  nearest  Cu-centered 
cluster.  From  crystallography,  the  center  of  a  tetrahedral 
site  in  an  FCC  lattice  is  located  at  a  distance  rT1  =  0.433ac 
of  the  nearest  cell  atom  (Cu-centered  cluster  in  our  case) 


r(A) 


FIG.  14.  Correlation  of  the  Cu  Cu  distances  in  the  FCC  super  lattice  formed 
of  the  Cu  centered  clusters  and  the  intensity  peaks  on  the  Cu  Cu  partial  ra 
dial  distribution  function. 


TABLE  XIV.  The  lattice  parameter,  ac,  the  number  of  atoms  per  unit  cell, 
Nc,  and  the  number  of  occupied  interstitial  sites  per  unit  cell,  V,  in  the  face 
centered  cubic  lattice  formed  of  the  Cu  centered  clusters  in  Ca  Mg  Cu 
amorphous  alloys. 


Alloy 

Ca4oMg25Cu35 

Ca45Mg25Cu3o 

Ca5oMg25Cll25 

Ca60Mg25Cu15 

ac  (A) 

10.9  ±  0.1 

1 1.0  ±  0.1 

10.9  ±0.1 

11.2  ±0.2 

Nc 

51.6  ±  1.5 

49.6  ±  1.5 

45.3  ±  1.3 

43.6  ±2.3 

N, 

6.4  ±  1.5 

4.0  ±  1.5 

0.0  ±1.3 

1.2  ±2.3 

and  rT2  =  0.83ac  of  the  second  nearest  cell  atom,  while  the 
center  of  the  octahedral  site  is  located  at  the  distance 
ro  =  0.5ac  of  the  nearest  cell  atom.  The  two  characteristic 
distances  between  the  octahedral  and  tetrahedral  sites  are 
ron  =  0.433ac  and  rOT2  =  0.83ac.  Therefore,  at  ac=  11.0  A, 
rT/  =  ron  =  4.8  A,  and  /'0  =  5.5A  correspond  to  the  second 
PRDF  peak,  and  rT2  =  rOT2  =  9.13  A  likely  corresponds  to 
the  fourth  Cu-Cu  PRDF  peak.  Figure  14  shows  good  corre¬ 
lation  between  the  Cu-Cu  distances  in  the  FCC  super¬ 
lattice,  formed  of  the  Cu-centered  clusters,  and  the  peak 
maxima  on  the  Cu-Cu  partial  radial  distribution  function.  It 
can  also  be  seen  from  this  figure  that  the  Cu-Cu  MRO  cor¬ 
relations  are  only  present  within  one  FCC  super-cell  and 
they  disappear  at  distances  higher  than  ac.  Similar  analysis 
was  conducted  for  the  other  three  Ca-Mg-Cu  glasses  and 
Table  XIV  provides  estimated  values  of  the  FCC  lattice  pa¬ 
rameter,  ac,  as  well  as  the  number  of  atoms,  Nc,  and  the 
number  of  occupied  interstitial  sites,  /V,,  per  unit  super-cell 
in  the  studied  amorphous  alloys. 

The  interstitial  sites  can  tentatively  be  occupied  not  only 
by  Cu  but  also  Mg  and/or  Ca.  Indeed,  a  detailed  analysis, 
similar  to  the  presented  above,  has  led  to  a  conclusion  that 
the  second  peak  in  the  Mg-Cu  PRDF  can  be  explained  by  the 
presence  of  Mg  inside  the  tetrahedral  and  octahedral  sites. 
However,  the  second  peak  in  the  Ca-Cu  PRDF  can  be 
explained  by  the  presence  of  Ca  inside  octahedral  sites  only. 
There  are  no  peaks  in  the  Ca-Cu  PRDF  corresponding  to  the 
tetrahedral  positions.  This  observation  can  be  easily 
explained  by  a  large  size  of  Ca  atoms,  which  would  not  fit 
the  small  tetrahedral  sites.66 

F.  Topological  description  of  Ca-Mg-Cu  atomic 
structures 

The  present  results  give  the  most  detailed  experimental 
picture  of  a  ternary  metallic  glass  to  date,  and  so  it  is  worth 
analyzing  the  current  findings  with  the  efficient  cluster  pack¬ 
ing  (ECP)  model  for  the  atomic  structure  of  metallic 
glasses.42  43  67  In  binary  glasses,  preferred  radius  ratios  and 
strong  chemical  interactions  between  unlike  atoms  produce 
efficiently  packed  solute-centered  clusters  that  are  arranged 
in  space  with  an  FCC-like  structure  over  a  length  scale  of  a 
few  cluster  diameters.  In  forming  these  clusters,  solutes  first 
tie  up  the  other  atoms  in  the  structure,  and  higher  solute  con¬ 
centrations  give  structures  where  the  cluster-octahedral  and 
cluster-tetrahedral  interstices  are  progressively  filled.  At 
even  higher  solute  concentrations,  solutes  replace  solvent 
atoms  in  the  cluster  1st  shell.  Structures  where  all  of  the 
cluster-interstitial  sites  and  about  1/3  of  the  1st  shell  are 
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occupied  by  solvent  atoms  generally  give  the  most  stable  bi¬ 
nary  glasses.  A  range  of  experimental  and  computational 
data  supports  this  model. 

It  was  earlier  suggested  that  the  largest  solute  in  a  multi- 
component  glass  would  be  the  primary  structure-forming  sol¬ 
ute,  since  larger  solutes  are  topologically  more  potent.67 
However,  the  present  work  supports  and  earlier  sugges¬ 
tion16'63  that  elements  with  the  strongest  chemical  interaction 
with  solvent  atoms  will  form  the  structure-forming  clusters. 

While  much  of  the  earlier  analyses  depend  upon  the 
nominal  radius  ratio  between  atoms,  it  has  more  recently 
been  shown  that  the  weighted  average  radius  of  atoms  in  the 
first  coordination  shell  gives  a  good  indication  of  GFA.68 
This  effective  radius  ratio,  r,  correlates  strongly  with  good 
GFA  and  gives  an  equivalent  structural  description  regarding 
local  partial  coordination  numbers.66  This  concept  is  used 
here  to  successfully  model  triplet  angles.  Although  the  nomi¬ 
nal  radius  ratio  for  Cu  and  Ca  is  about  0.63  supporting  a 
CN  =  8,  the  effective  radius  ratio  is  about  0.72.  This  suggests 
efficiently  packed  clusters  with  CN  =  9,  which  is  less  than 
the  value  of  ~  10  in  the  present  work.  The  cut-off  distance  in 
the  present  work  is  about  30%  longer  than  rmode  for  Ca  and 
Cu  and  50%  larger  for  Mg,  which  may  contribute  to  this 
discrepancy. 

The  present  data  show  different  site  filling  rules  for  ternary 
glasses.  Inter-cluster  sites  are  essentially  vacant  in  Ca-Mg-Cu 
glasses  with  Cu  concentrations  of  15%  and  25%  (Table  XIV), 
but  Cu-Cu  contact  occurs  in  all  glasses  (Table  II),  suggesting 
that  some  Cu  atoms  occupy  sites  in  the  1st  coordination  shell 
before  inter-cluster  sites  are  filled.  Topological  analysis  sug¬ 
gests  that  sufficient  Cu  exists  to  begin  filling  some  cluster- 
interstitial  sites  at  15%  Cu,  and  the  absence  of  this  structural 
feature  is  consistent  with  the  presence  of  Cu  atoms  in  the  1st 
shell  of  Cu-centered  clusters. 

Solute-rich  binary  metallic  glasses  with  solute  atoms  in 
the  1st  coordination  shell  of  solute-centered  clusters  are 
expected  to  have  4  cluster-interstitial  atoms  per  structure¬ 
forming  cluster,42’43’67  or  16  cluster-interstitial  atoms  per 
FCC  unit  cell.  However,  recent  phenomenological  analyses 
of  multi-component  BMGs  show  that  they  contain  only  1  or 
3  solute  atoms,  or  “glue  atoms,”  per  cluster.68’69  In  agree¬ 
ment  with  this,  the  present  work  shows  1  cluster-interstitial 
site  per  Cu-centered  cluster  in  Ca45Mg25Cu3o  and  1.5  sites 
per  cluster  for  Ca40Mg25Cu35. 

With  this  as  background,  the  present  data  give  the  follow¬ 
ing  topological  picture  of  Ca-Mg-Cu  structures.  Cu-centered 
clusters  in  Ca6oMg25Cu15  have  a  total  CN  of  9.8,  so  that  the 
bulk  concentration  of  Cu  is  achieved  with  1.7  Cu  atoms  per 
cluster.  There  are  no  Cu  atoms  in  cluster-interstitial  sites 
(Table  XIV),  so  there  is  1  Cu  atom  at  the  center  of  the  cluster 
and  0.7  Cu  atoms  per  cluster  in  the  1st  coordination  shell.  This 
agrees  very  well  with  CNCu_Cu  of  0.6  (Table  II).  The  composi¬ 
tion  of  Cu-centered  clusters  in  Ca60Mg25Cu15  is  nearly  equal 
to  the  bulk  composition,  further  supporting  this  interpretation. 
The  total  CN  around  Cu  clusters  is  10.4  in  Ca50Mg25Cu25,  so 
there  are  3.5  Cu  atoms  per  structure-forming  cluster.  Data  in 
Table  XIV  suggest  there  are  no  atoms  in  cluster-interstitial 
sites,  so  that  there  would  be  1  Cu  atom  in  the  cluster  center 
and  2.5  Cu  atoms  in  the  1st  coordination  shell.  However,  this 


gives  a  partial  CN  that  is  larger  than  the  value  in  Table  II,  and 
the  composition  of  Cu-centered  clusters  is  noticeably  deficient 
in  Cu.  Since  the  data  in  Table  XIV  have  the  largest  uncer¬ 
tainty,  this  structural  interpretation  is  unlikely.  The  bulk  com¬ 
position  is  achieved  when  0.7  Cu  atoms  in  cluster-interstitial 
sites  are  added  to  the  Cu-centered  cluster  composition,  so  that 
an  alternate  interpretation  is  that  there  is  1  Cu  atom  at  the  cen¬ 
ter  of  the  clusters,  0.7  Cu  atoms  in  cluster- interstitial  sites,  and 

1.8  Cu  atoms  in  the  1st  coordination  shell  of  the  Cu-centered 
clusters.  This  structure  gives  a  partial  CNCu_Cu  that  is  consist¬ 
ent  with  the  data  in  Table  II  and  with  the  observation  of  a  peak 
in  the  Cu-Cu  g(r)  near  7.8  A.  The  total  Cu  CN  is  10.4  in 
Ca45Mg25Cu3o,  so  there  are  4.5  Cu  atoms  per  cluster.  The 
composition  of  Cu-centered  clusters  is  deficient  in  Cu,  so  that 
an  additional  0.7  Cu  atoms  per  cluster  are  needed  to  match  the 
bulk  composition.  This  suggests  1  Cu  atom  at  the  cluster  cen¬ 
ter,  0.7  Cu  atoms  at  cluster-interstitial  sites  (Table  XIV),  and 

2.8  Cu  atoms  in  the  1st  coordination  shell.  This  is  consistent 
with  the  data  in  Table  XIV,  with  bulk  composition  and  partial 
CNs,  and  with  the  Cu-Cu  peak  in  g(r)  near  7.8  A.  Finally,  there 
are  5.5  Cu  atoms  per  structure-forming  cluster  in 
Ca40Mg25Cu35.  To  match  bulk  composition,  1.8  Cu  atoms 
must  occupy  cluster-interstitial  sites.  Thus,  there  is  1  Cu  atom 
at  the  cluster  centers,  1.8  Cu  atoms  in  cluster-interstitial  sites, 
and  2.7  Cu  atoms  in  the  1st  coordination  shell.  This  structure 
is  consistent  with  the  number  of  occupied  interstitial  sites  from 
Table  XIV  (1.6  ±  0.3),  with  bulk  composition  and  Cu-Cu  par¬ 
tial  CNs  (Table  II),  and  with  the  Cu-Cu  g(r). 

While  bulk  compositions  can  be  matched  by  adding  Cu 
atoms  to  cluster-interstitial  sites,  the  Mg  concentrations  are  still 
slightly  deficient.  The  present  work  suggests  that  some  Mg  is  in 
cluster-interstitial  sites,  consistent  with  this  observation. 

V.  CONCLUSIONS 

(1)  The  atomic  structures  of  four  C a4 0 + x M g 2 5 C u  3  5  x 
(X  =  0,  5,  10,  and  20  at.  %)  ternary  metallic  glasses  have 
been  determined  using  a  synergistic  combination  of  neu¬ 
tron  diffraction,  ab  initio  molecular  dymamics  (QMD) 
and  CRMC  modeling.  All  six  PRDFs,  gij(r),  have  been 
identified  for  each  alloy.  The  S(Q)  curves  damp  very 
quickly  with  Q  and  no  oscillations  are  apparent  beyond 
~15A  1 ,  indicating  that  the  glasses  have  a  wide  range 
of  interatomic  distances.  An  increase  in  Cu  content 
shortens  the  average  interatomic  distance  and  narrows 
the  distribution  of  interatomic  distances,  shifting  the  first 
sharp  diffraction  peak  to  higher  Q  and  broadening  the 
peak.  A  pre-peak  at  Q«  1.2  A  1  in  the  S(Q)  curves  indi¬ 
cates  the  presence  of  MRO  in  these  amorphous  alloys. 

(2)  The  first  PRDF  peaks,  which  correspond  to  interatomic 
bond  distances  in  the  1st  coordination  shell,  are  asym¬ 
metric,  and  the  mode  bond  distances  are  always  smaller 
than  the  respective  mean  bond  distances.  The  difference 
is  the  smallest  for  Ca-Ca  pairs  and  increases  in  the  order 
of  Ca-Mg,  Mg-Mg,  Ca-Cu,  Cu-Cu,  and  Mg-Cu. 

(3)  The  nearest-neighbor  mode  bond  lengths  are  shorter  than 
those  in  competing  crystals.  A  noticeable  shortening  of 
Ca-Cu  and  Mg-Cu  bond  distances  indicate  strong  inter¬ 
actions  between  these  atom  pairs.  It  is  suggested  that  the 
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bond  shortening  is  enabled  by  the  absence  of  long-range 
atomic  order,  which  lowers  the  free  energy  of  metallic 
glasses  and  increases  GFA. 

(4)  Pronounced  CSRO  near  Cu  atoms,  CSRC  near  Ca  atoms, 
and  a  neutral  environment  near  Mg  atoms  are  shown. 
Increasing  the  Cu  concentration  from  15  to  35  atomic 
percent  increases  the  total  coordination  number  around 
Ca,  CNCa,  from  13.6  to  15.0,  while  CNMg  and  CNCu 
remain  unchanged  at  those  «12.5  and  «  1 0.4.  respec¬ 
tively.  The  partial  coordination  numbers  depend  nearly 
linearly  on  the  concentration  of  Cu,  so  that  the  number 
of  Ca  atoms  decrease,  the  number  of  Cu  atoms  increase 
and  the  number  of  Mg  atoms  are  almost  constant  with 
increasing  Cu  concentration. 

(5)  Voronoi  tessellation  and  three-body  correlations  show  that 
many  types  of  coordination  polyhedra  are  present,  but  the 
most  common  are  (0,2, 8, 4)  Ca-centered,  (0,2, 8, 2)  Mg- 
centered,  and  (0, 3,6,0)  and  (0,2,8, 1)  Cu-centered.  The 
fractions  and  distributions  of  these  clusters  depend  on 
alloy  composition.  Polytetrahedral-type  clusters  and  live- 
coordinated  vertices  dominate  in  the  amorphous  struc¬ 
tures,  which  indicate  that  tetrahedra  and  pentagonal 
bi-pyramids  are  the  main  building  blocks  in  these  amor¬ 
phous  alloys.  The  GFA  of  the  Ca-Mg-Cu  alloys  increases 
with  a  decreasing  fraction  of  Cu-centered  (0,3,6,0)  clusters. 
The  local  atom  packing  fraction  near  a  given  atom  increases 
with  atomic  radius  and  with  an  increase  in  Cu  concentra¬ 
tion.  Global  atom  packing  fractions  of  0.68  to  0.73  are 
much  larger  than  for  dense  random  packed  structures  and 
are  comparable  to  packing  in  crystalline  structures. 

(6)  Medium  range  order  is  present  in  Ca-Mg-Cu  amorphous 
alloys  and  is  the  result  of  strong  tendency  of  Cu  for 
CSRO.  Detailed  analysis  shows  that  the  amorphous 
structure  of  these  alloys  is  described  by  close-packing  of 
Cu-centered  clusters  that  follows  a  face-centered  cubic 
motif  over  a  length  scale  of  ^10  A  or  about  1.5  cluster 
diameters.  This  gives  a  characteristic  medium  range 
order  pre-peak  at  Q~1.2A  1  in  the  total  scattering 
structure  factors  of  these  alloys. 

(7)  A  topological  description  of  Ca-Mg-Cu  atomic  structures 
shows  that  the  effective  radius  of  atoms  in  the  1st 
coordination  shell  of  Cu  atoms  enables  efficient  atomic 
packing.  These  clusters  match  the  bulk  composition  for 
CaA(lMg25Cu  1 5  but  are  Cu-lean  for  all  other  compositions. 
The  remaining  structures  have  between  0.7  and  1.8  Cu 
atoms  in  cluster-interstitial  sites.  This  allows  the  bulk  com¬ 
position  to  be  matched  and  is  necessary  for  agreement 
with  the  Cu-Cu  g(r)  plots.  Cu  atoms  occupy  between  0.7 
(for  Ca60Mg25Cu15)  and  2.8  (for  Ca45Mg2gCu3o)  sites  in 
the  first  coordination  shells  of  Cu-centered  clusters. 
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